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1. INTRODUCTION 


One of the major difficulties associated with the 
design and use of composite materials for load-bearing struc- 
tural applications is the current lack of a suitable failure 
criterion for both the individual laminae and the laminated 
structure as a whole. Although many lamina failure criteria 
have been proposed {see Refs. 1, 2 and 3, for example), in- 
sufficient experimental data particularly under combined 
states of stress have been accumulated to -indicate which cri- 
terion is best able to predict the failure stresses. The 
difficulty is of course that each strength criterion has been 
developed empirically with certain interaction parameters 
being neglected and they are all phenomenological in nature. 
In other words, they can predict the occurrence of failure 
but they do not describe the physics or mode of failure. 

It would appear that the most general failure cri- 


terion proposed up to the present is that given by Wu {Refs. 
4, 5, 6) in the form of a tensor polynomial 


f^i) 


f. cr. + f. .* cr. or. 

ii id i D 


O'. <r. cr. 

i D k 


t .... — 1 


( 1 ) 


(where i, j, k = 1, 2, 3, ...6) 

which can be shown to encompass all other failure criteria 
which are currently available. The simplest form of Eq. (1) 
which retains the interaction tensor strength components is 
(Refs. 7, 8) 
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( 2 ) 
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This quadratic tensor polynomial defines a failure surface 
in stress space in terms of two strength tensors and 
of the second and fourth ranks, respectively. Of particular 
interest in this formulation is the existence of the linear 
terms in k) which can account for the observed differences 
between positive and negative stress induced failures. How- 
ever, in order to employ the stress tensor polynomial strength 
equation, one is faced with the difficulty of evaluating not 
only the standard principal strength parameters as defined 
by the (F^, F^) relations, but also determining the inter- 
action terms , F ijk' etc ') which are regarded as inde- 

pendent material properties. Since the failure surface may 
not be ellipsoidal in shape {i.e., the principal directions 
of strength may not always be orthogonal) , it is necessary 
to include higher order terms in the tensor polynomial equa- 
tion (such as the sixth-order failure tensor F. .p) . Thus the 
number of independent strength parameters that have to be de- 
termined experimentally can become inordinately large. It 
might be noted that a compromise "hybrid" method has ( been re- 
cently proposed (Ref. 6) in which fewer experimentally mea- 
sured interaction parameters are required. One other formu- 
lation that has also been published recently (Ref. 9) is based 
on strain energy and matching the shape of the failure surface 
to experimental data. However, only the three principal 
strength parameters are utilized and no interaction terms are 
included. 
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Most of the experimental studies that have been under- 
taken with composite materials have employed flat plate speci- 
ments , particularly for evaluating the orthotropic stiffness 
and principal strength parameters. For example, two reports 
by Azzi and Tsae (Refs. 10, 11) have demonstrated remarkably 
well how the plane stress anisotropic Hill failure condition 
was able to predict the strength for transversely isotropic 
configurations, including unidirectional off-axis, cross and' 
angle (+0) ply laminates. In this failure criterion, only 
the primary strength parameters were used, with no allowance 
for compression/tension differences. Since theory and experi- 
ment were in good agreement, one must conclude that for the 
specimens tested, the interaction terms must either have 
been reasonably well approximated or their influence on the 
strength prediction was of second order. However, in general 
this may not be the case for all materials and load config- 
urations. In particular, for the flat off-axis test speci- 
men, Pipes and Cole (Ref. 12) recently showed very well that 
the F ^2 term did not contribute significantly to the tensile 
strength prediction as a function of fiber orientation (0). 
Moreover, it was proven that this specific test configuration 
was not sufficiently sensitive to accurately measure F^ 2 . 

Other difficulties associated with flat plate specimens include 
- specimen geometry can affect the measured strength - thus 
it is necessary to use samples having very high (free length/ 
width) ratios (Ref. 13)' 
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existence of free edges can lead to large interlaminar 
shear strains and premature cracking of matrix (Ref. 14); 

- effect of end constraints and coupled deformations can 
lead to erroneous results; 

obtaining a homogeneous stress state requires a sufficiently 
wide specimen, but this makes it difficult to overcome the 
first problem. 

Consequently, a great deal of theoretical and experi- 
mental work has been devoted to assessing the circular cyl- 
indrical tube as an optimum specimen configuration for both 
stiffness and strength characterization. However, consider- 
able care must also be exercised in the proper design of the 
cylinder to ensure absence of edge constraint effects, fail- 
ure in the test section and a uniform stress field. In gen- 
eral, to employ "thin wall" shell theory for calculating 
stress levels requires that the (diameter/thickness) ratio 
should be > 20 (Refs. 15, 16) and the specimen length should 
equal twice the mean diameter plus the desired gauge length 
• 15) . The limits on (-diameter/thickness) however are 
dependent on the material system and ply orientations. It 
was shown, for example, that it was difficult to achieve a 
uniform stress distribution in a helical unidirectional wound 
tube and thus one should use + 0 symmetric configurations 
fo-c characterization tests. Perhaps the most difficult para- 
meter to assess experimentally is the degree of end attach- 
ment needed to obtain a relatively uniform stress distribution 
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in the gauge section without failure occurring at the ends 
(Refs. 17, 18) . For a good review of current static test- 
ing techniques, one can refer to Bert's paper (Ref. 19). 

At this point, it is worthwhile reviewing the extent 
to which tubular specimens have been utilized to evaluate 
strength criteria. Aside from the preliminary burst press- 
ure results published by Sandhu et al (Ref. 20), the first 
measurements of an interaction component (F.^) were made by 
Wu (Ref. 8). More details on the actual test programme can 
also be found in Ref. 21. Later, when Wu expanded his fail- 
ure criterion to include higher order terms (F^) / he also 
reported additional data (Ref. 6) in the form of a failure 
surface for direct comparison with various failure theories. 
These results were also given in Ref. 5 in which it was shown 
that the tensor failure polynomial provided the best compari- 
son. However, the scatter was substantial and the configura- 
tions tested were quite restrictive, i.e., specially ortho- 
tropic and symmetric. One noteworthy investigation involved 
the testing of + © cylinders under axial tension, torsion and 
internal pressure (Ref. 22) . The failure criterion adopted 
consisted of evaluating the in-plane stresses and determining 
when their fracture stress levels had been reached (based on 
four fracture stress parameters - tensile and compressive 
stresses for fibers, tensile stress normal to fibers and 
shear stress) . An empirical expression was used to calculate 
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the ultimate shear stress as a function of orientation with 
respect to the fibers. Furthermore, it was assumed that the 
matrix failed first and the stresses could be re-calculated 
based on an assumed zero modulus for the matrix. With this 
approach, good agreement between experiment and predicted 
behaviour was achieved. Whether or not this method would 
yield similar results for other laminates is not known. Fin- 
ally, one last report which was recently published (Ref. 23) 
contains experimental data on tubes subjected to uniaxial 
loading, internal pressure and torsion. Some comparisons 
with failure theories were presented but no attempt to eval- 
uate the tensor strength polynomial terms or its strength pre 
dictions was made. 

The following report describes the experimental mea- 
surements and techniques used to obtain the strength tensor 
components, including cubic terms. Based on a considerable 
number of biaxial pressure tests together with specimens sub- 
jected to a constant torque and internal pressure, a modi- 
fied form of the plane stress tensor polynomial failure equa- 
tion was obtained that was capable of predicting ultimate 
strength results reasonably well. Although only glass-epoxy 
tubular test configurations have been studied, it is felt 
that the failure' criterion should apply equally well to other 
material laminates. In addition, preliminary data were ob- 
tained to determine the effect of varying post cure times 
and ambient temperature (- 80° F to 250° F) on the change in 



7 


two tensor strength terms, F 2 and F 22 . Other laminate con- 
figurations presumably should yield corresponding variations 
for the remaining strength parameters. 
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2 . FABRICATION OF SPECIMENS 

2 .1 Belt Wrapper ■ 

The belt wrapper, as shown in Figure 1, is an apparatus 
which applies constant pressure around approximately 340° of 
a mandrel through the use of a silicone coated fibreglass 
belt under tension. The mandrel is positioned on the loose 
silicone/fibreglass belt between two rollers. One of these 
rollers can be moved forward in a groove and tightened down 
so that there is a gap of about 1,5 inches between the two 
rollers (Figure 2), One of the lower rollers is mounted at 
both ends on air cylinders and these cylinders are connected 
through a pressure regulator to an air supply. When the air 
cylinders are pressurized, the lower roller is pushed for- 
ward so that tension is applied to the belt. This causes 
the mandrel to be pressed up against the two upper rollers 
and the belt applies pressure around the mandrel. There is 
a reversible, variable speed motor connected to one of the 
lower rollers. When the motor is engaged, the roller rotates 
causing the belt to move, which, in turn, causes the mandrel 
to rotate. In this way, the prepreg tape can be wrapped 
tightly onto the mandrel. 

The mandrels can be made from either steel or aluminum 
tubing of the desired outside diameter. The tube must be 
cut to a length of 46 inches so that it will fit in this 
particular belt wrapping apparatus and then, its surface must 
be ground to ensure a uniform circular cross-section along 
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the entire length and be. very smooth. The tube is 
then fitted with a vacuum connection at one end. Now, the 
finished mandrel is thoroughly cleaned to remove all traces 
of dirt and oil from its surface and an even continuous film 
of Frekote 33 Releasing Interface is sprayed on. This film 
is then baked on at 125 *C for about 15 minutes to increase 
its durability as suggested by the manufacturer. 

2,. 2 Geometry of Pre-Preg Sheets for Wrapping 

In order to wrap a particular fibre angle on the mandrel 
and to avoid overlapping or gaps at the seams of the material, 
the pre- impregnated (pre-preg) tape must be cut very accu- 
rately to the desired width. The following procedure is used 
to determine the length and width of the tape required to 
fabricate a specimen having any fibre angle, © . 

Figure 3 shows the dimensions of the finished specimen 
and of the tape required to make that specimen. In this 
figure, D is the outside diameter of the mandrel, L is the 
required length of the test specimen, C is the outside cir- 
cumference of the mandrel, W is the width of the pre-preg 
tape, 1 is the length of the pre-preg tape and © is the fibre 
angle. The dimensions of the pre-preg tape are calculated 
using the following relations: 

C = ttD, W = ttD sin ^ , and 1 = L / cos © 

where ^ = 90 - © degrees. When performing these calculations 
for additional laminae other than the first, the value of D 
must be increased by twice the thickness of the pre-preg 
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tape each time. 

2.3 Wrapping; Procedure 

The material used to make all of the test specimens was 
"Scotchply" Reinforced Plastic Type 1002 f ibreglass/epoxy 
pre impregnated tape. 

The feed guides on the belt wrapper are adjusted to 
produce butt seams using a paper pattern of the desired ori- 
entation. The first ply of pre-preg is then placed on the 
working table of the belt wrapper with the paper backing 
against the belt. The circumferential edge of the tape is 
positioned against the feed guide and the material is slowly 
rolled onto the mandrel in such a way that the paper backing 
faces outwards. If the tape has been cut correctly and the 
guides adjusted carefully, there should be no overlaps or gaps 
in the tube, but rather a smoothly butted seam. The paper is 
then peeled away from this first ply by running the belt 
wrapper backwards and carefully lifting the paper from the 
tape. Additional laminae are wrapped in the same manner. 

With thick-walled structures, the air cylinder pressure should 
be reduced as successive layers are wrapped to avoid buckling 
the material. 

Once all of the plies are wrapped onto the mandrel, a 
porous teflon coated fibreglass cloth is wrapped around the 
tube again using the belt wrapper to ensure a wrinkle-free 
application. The specimen is then removed from the belt- • 
wrapper and bagged in Vac-Pak type E3760 film and sealed 
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with a vacuum bag sealant# Canvas strips are placed along 
the bag seam and around the ends of the mandrel, inside the 
vacuum bag, to allow a vacuum over the entire tube (Figure 4) 
Using a vacuum pump, a vacuum of 29 inches of Hg is es- 
tablished and maintained in the bag around the specimen for 
three hours while the tube is curing at 150 W C. The bag and 
release cloth are removed from the tube after four hours and 
the specimen is post cured for seventeen hours at 150®C. The 
oven temperature is then raised to 200°C for forty-five 
minutes* After this, the tube is removed from the oven and 
slid off the mandrel. 

The tubes are cut to the desired length by sliding them 
over a cutting mandrel mounted on a lathe (Figure 5). Using 
a slow lathe speed, the tube is turned and an air operated 
abrasive cutting disc, mounted in the tool post, is used to 
make a square cut. The tube ends are sanded smooth and thick 
ness .measurements are taken at eight positions equally spaced 
around the circumference at both ends and in the middle 
section in order to obtain an average tube wall thickness. 

The specimens are then readied for testing by reinforcing 
the ends with stepped down layers of cotton cloth, fine fibre 
glass cloth and epoxy and potting them in end plates, made 
of a suitable material, to a depth of one inch with a room 
temperature curing epoxy. Several specimens were strain 
gauged with 350 ohm, 0.5 inch foil gauges to provide stress- 
strain data up to failure. 
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2 . 4 Specimen Design 

Thin-walled cylindrical test specimens offer definite 
advantages for experimentally characterizing the mechanical 
behaviour of fibre/matrix composite materials. All of the 
loadings required to fully characterize the composite system 
with respect to its mechanical properties can be applied to 
this type of specimen. However, particular care must be taken 
to ensure that reasonably uniform stresses across the wall 
thickness are obtained. As noted earlier, the stresses in- 
duced in an anisotropic tube which is subjected to combined 
axial load, torsion, and internal pressure are approximately • 
uniform across the wall thickness provided that the ratio of 
wall thickness to radius is sufficiently small. 

A major problem remains in the use of tubular test 
specimens which are usually rigidly clamped {i.e., potted) 
to end plates. This method of end attachment prohibits the 
transverse and radial displacements of the tubular specimen 
and this induces relatively large bending stresses in the 
region near the end plates. - These bending stresses can lead 
to premature failure of the tube in the end regions. 

If the magnitude of these bending stresses can be pre- 
dicted before any testing is done, then sufficient reinforce- 
ment can be applied to the ends of the tube in order that 
they do not fail prematurely. Also, if the variations in the 
induced stresses along the length of the test specimen are 
known, then it can be designed more accurately so as to 
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ensure a uniform stress distribution in the test section. 

i 

In Appendix A, the governing equations for ortho- 
tropic symmetric laminates (A lg = A 2g = B^.. = 0) under symm- 
etric loading are derived using an extension of Flugge's 
shell theory (Ref. 24) provided by Cheng and Ho (Ref. 25). 
These equations were then solved to determine the proper dim 
ensions of the tubular test specimens used for composite 
characterization . 
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3. EVALUATION OF STRENGTH TENSORS 

The general form of the tensor polynomial strength 
equation proposed by Wu (Ref. 4) is 

F. <T . + F..CT. T . + F . .y, . <y. CT- + ••• =1 (3) 

where i, j, k = 1, 2, ... 6, and F^, F^, are tensors 

of the 2nd, 4th and 6th ranks, respectively. For a plane 
stress state, i, j, k = 1, 2, 6. The purpose of this section 
of the report is to provide a description of the experimen- 
tal methods used to evaluate the principal strength compon- 
ents (F^, , Fg, F.^, F 22' F 66^ anc ^ < 3 uac ^ rat i c interac- 
tion parameters F -^g a nd F^g) , assuming the material 

system is symmetric (i.e., F^ ^ = F^ for i f j). In this 
formulation, each of the tensor strength terms (F^, F ^j/ 

F. .p, etc.) is considered to be an independent material pro- 
perty. Although experimental results have only been obtained 
for the linear and quadratic terms, analysis in Section 6 was 
used to estimate the required cubic terms (F. ) necessary 

for the strength equation. 

3.1 Principal Strength Tensor Components (F^, F^) 

From the analysis by Wu (Ref. 8), it was shown that 
the principal strength -tensor components (F^ and F^) can be 
readily calculated from the experimentally determined values 
of the uniaxial tensile and compressive failure stresses in 
the fiber direction (X and X 1 ) , perpendicular to the fibers 
(Y and Y ( )and from positive and negative pure shear failure 
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stresses (S and S', respectively). The appropriate relations 
are given by: 



F=-_± p - 

6 S S' 'll XX' (4) 

F 22 = W F 66 “ SS' 

where the terms are not considered since they were shown 

in Ref. 6 to be redundant and not necessary for the strength 
criterion. 

To determine the tensile strength in the fibre direc- 
tion, four tubes were tested. These tubes all had an inside 
diameter of approximately 2 inches and were 8 inches long be- 
tween end-plates. They were fabricated from two plies of 
material with the fibers aligned along the longitudinal axis 
of the tube, (i.e., © =0). Figures 6 to 8 show the test 

set-up used and a specimen before, during and after failure. 
The results are summarized in Table I. A typical stress-' 
strain curve, obtained using a strain gauge oriented longi- 
tudinally at the midpoint of the specimen, is shown in Figure 
9. 

For the compressive strength in the fiber direction, 
seven tubes were tested. It was necessary to fabricate these 
tubes with a smaller inside diameter (approximately one inch) 
and a larger wall thickness in order to avoid buckling. They 
were made from eight plies of material and the fibers were 
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also aligned in the longitudinal direction. The test set- 
up and a failed specimen are pictured in Figures 10 and 11. 
Figure 9 also shows a stress-strain curve for one of the 
tubes . 

Seven unidirectional, four ply, laminated tubes with 
& = 90 ° were tested in tension to determine the value of the 
ultimate transverse tensile strength. These tubes had an 
inside diameter of about 2 inches and were 5 inches long be- 
tween end-plates. They were wrapped from a single piece of 
the pre-preg tape in order to provide continuous fibres in 
the circumferential direction. This resulted in two seams 
opposite each other on the inside and the outside walls of 
the tube. Although a slight discontinuity occurs at the seam, 
it does not affect the determination of the transverse 
strength at all. A tube is shown in the test set-up, after 
tensile failure, in Figure 12. The tensile strength results 
are contained in Table I and Figure 14 shows the stress- 
strain behaviour of one of the test specimens. 

The compressive strength of the material transverse 
to the fiber direction, was determined by applying axial com- 
pression to failure on seven specimens, similar to those used 
in the transverse tension tests. Figure 13 shows a typical 
failure of one of the tubes. These results are also summar- 
ized in Table I and Figure 14 shows a stress-strain curve 
for this type of loading. 


/ 



17 


Seven tubes were tested to failure under pure tor- 
sional loading to obtain the ultimate shear strength of the 
material. As expected, the direction of the torque did not 
significantly change the failure stress. The tubes were 
four-ply laminates with 6 = 90°, fabricated in the same way 
as those that were used in the tests for the transverse ten- 
sile and compressive strengths. Strain gauges were applied 
at 0 = + 45° to some of the tubes so that the stress-strain 
data to failure could be obtained. One of these curves is 
shown in Figure 15. The shear failure data are presented in 
Table I and Figures 16 and 17 show the test set-up and a 
failed tube, respectively. 

Once the five basic failure strengths of the material 
are known, the simple strength tensors can be easily calcu- 
lated using the relations given in Eq. (4). These values 

are listed in Table II. Using the stability condition (Ref. 

2 

7) F jj “ F ^j > 0, the bounds on the interaction strength 

tensors can also be calculated (Table II) . 


3.2 Quadratic Interaction Strength Tensor Components (F,^..) 
It was shown in the previous section that the deter- 
mination of the F^ and F^ (i = 1, 2, 6) strength tensors re- 
quires only simple uniaxial tension, compression and torsion 
tests. On the other hand , the interaction tensors F._, F. , 

^ X D 

and require biaxial stress experiments in which very 

strict control must be maintained over the biaxial stress 
ratio, B = CT. /0T. in order to achieve sufficient test accuracy 
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to correctly measure the F ± . terms. In Ref. 8, Wu has shown 
that the best resolution of the interaction tensors is ob- 
tained by using an optimal value of the biaxial stress ratio 
together with the most suitable stress state since the re- 
solution of the interaction tensor is dependent on both of 
these variables. However, the major difficulty that imme- 
diately arises is the fact that the full strength equation (1) 
must be employed, thus introducing cubic (and higher order) 
terms. As noted earlier, we shall restrict our attention 
only to a cubic formulation. 

Consider, for example, the relative scatter (or resolu- 
tion) of the interaction tensor component From Ref. (6) 

this is given by 


A F tt ^ 




A<^ 


(5) 


where 0” 2 is the failure stress of a biaxial test at a given 
stress ratio of B 12 = ^ ^ 2 is obtained as a toot of 

the failure equation (1) which can be re-written in the form 
(Ref. 6), 


( 3 F|,2 6,2 + 3 F|22 6,2 ) ^2 ~t ( Br + 2 
-h ^22 ^ ^2 2 +■ ( I &I2 ^ ^2 ) ^2 ~ ^ 


where CTg = 0 by definition. 

In order to alleviate the experimental requirements 
for determining F^ 2 , 


as a first approximation we shall assume 
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that the cubic contribution to the strength (i.e., ? 122 and 
^^^2 terms) is negligible for the 0° and 90° configuration. 
Consequently, one 'can then estimate F^ 2 and subsequently 
employ other complex load tests to assess the accuracy of 
its value. Based on these results, it is possible to estimate 
values for F 122 , F 112 and then reiterate to correct F 12 . This 
procedure can be repeated to arrive at final values for F, 9 , 


F 112 and F 122* 


Hence, neglecting the cubic terms one obtains 


dF, 


IZ 


where 


Z ~— l 


= - ( F„ B .z + 2 F iz + Ik 'll. + ^)^z (7) 


B 


iZ 


2 


^ - -( + F z) i [ ( F| B l2 + F 2 ) z + 4 ( F„ B,2 


F2F )2 B )2 t 1 ^ 2 )] / 2 (F,, B t2 t 2F a B |X +* F 2Z ) 


( 8 ) 


The optimum biaxial stress ratio B^ 2 occurs when 


d V 


\z 


d B 


>2 
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b®\2 

A F a 


dBa 




Substituting into Eq. (9) yields 
^2(F (| 8,^F a ) 4- F, _ 

^ Q 'z (F u - F 2Z ) - F 2 


(9) 


Fa " \ ^2 

scatter magnification factor (Ref. 6). 


~ C ^ ^ ( Fjj B i2 + 2 F )2 
tF a )+ F, B a +■ F 2 ] 

B R (B a F 1+ F 2 ) (10) 
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where O'* can be determined from Eq. (8) . To obtain the 
best estimate for 2 t a set of experiments must be performed 
since the optimal value for B^ 2 depends on F.^/ assuming the 
other tensor strength components are known. 

Before any strength tests could be carried out to de- 
termine the interaction tensor components, an analysis had 
to be done to determine the optimal biaxial stress ratios for 
each stress state. This was accomplished once the bounds on 
each of the interaction tensors were known (Table II). 

Again, for example, consider the F^ 2 interaction tensor. 
Once the constants , F 2 , F.^, F 22 an< ^ t ^ e bounds on F^ are 
determined, it is a simple matter to substitute these into 
Eqs . (8) and (10) to estimate the optimal biaxial stress 

ratio for various values of F. 0 within the bounds of the 
stability condition. There are four roots of Eq. (8) for each 
value of f ^ 2' one ^ or eac ^ stress quadrant. . Then, the resolu- 
tion of F^ for the optimal ratios B^must be plotted to de- 
termine in which stress quadrant the testing should be done. 

The resolution of F^ 2 is given by the following relation; 

.= -(F.,B a + 2F, 2 + F M t f, \ 

(3 j2 2 q- 2 2 8 a <T 2 (11 ^ 

When these two graphs have been obtained, it remains to decide 
on the method of ' application of the biaxial stress state and 
then to begin the iterative procedure to determine the value 
°f f 22 * This method is equally applicable to the determina- 
tion of F^g and F 2 g. 
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interaction Tensor 

Figures 18 and 19 are graphs of the optimal biaxial 
stress ratio to be used in the determination of F^ an d the 
attainable resolution of F^ 2 for a given optimal B^ 2 . From 
these figures it is evident that for F^ 2 > 0, a longitudinal 

tension-transverse compression ( < 3T^> 0,CT 2 -< 0) experiment 
with the ratio B^ 2 in the range of - 14.4 to - 15.6 is optimal, 
and a" longitudinal compression-transverse compression ( CT^ ■< 

0, <T 2 < 0) experiment with the ratio B^ 2 in the range 13.2 
to 15.4 is optimal for F^ 2 0. The tension-compression ex- 

periment was chosen because it was the most convenient test 
for a tubular specimen. A combination of internal pressure 
and axial compressive loading was used to provide the required 
biaxial stress state. 

The specimens tested were unidirectional three ply- 
laminates with © = 90°. The tubes were continuously wrapped 
from a single piece of pre-preg tape and had a 1/4 inch over- 
lap at the seam for load transfer, A 1/16 inch polyurethane 
liner was spuncast inside each specimen to prevent any damage 
that might be caused by the oil used in the internal pressure 
apparatus . 

A schematic diagram and pictures of the actual test 
set-up are shown in Figures 20 and 21. An air-operated 
hydraulic oil pump was used to maintain the reservoir at a 
pressure of approximately 6‘, 00 0 psi. The specimens were 
pressurzed from this reservoir by adjusting a flow valve 
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while monitoring the internal pressure by means of a pressure 
transducer downstream from the valve. The specimens were 
mounted in a Tinius Olsen 60,000 lb. Universal Testing Ma- 
chine which provided the axial compressive load. This load 
was converted to a voltage by means of a potentiometer mounted 
in the testing machine. By knowing the pressure transducer's 
calibration in psi/volt and the testing machine's calibra- 
tion in lb/volt, the required pressure versus axial com- 
pressive load curve can be plotted for a particular biaxial 
stress ratio This curve was recorded on an X - Y plotter 

having internal pressure and axial compressive load as inputs. 
The valve was opened just enough to allow a slow, steady in- 
crease in pressure in the tube, while the amount of axial 
load applied by the testing machine was controlled manually 
to ensure that the loading followed the pre-calculated load 
curve up to failure. In this way, a constant biaxial stress 
ratio was maintained throughout the test. Using the values 
of the internal pressure and compressive load at failure, 
the failure stresses^ and 0" 2 were calculated and a value 
for the F 12 tensor was determined from the following equation: 

F 12 = 1 - (F ll Cr i 2 + F 22 !T 2 2 + F 1 <T 1 + P 2 Cr 2>] < 12 > 

This method provided a controlled, constant ratio bi- 
axial loading as is evident in Figure 22 which shows the actual 
loading curve in comparison to the calculated loading curve 
for a biaxial stress ratio of - 14.2. The results of the 
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iteration procedure to determine the value of F^ 2 are listed 
in Table III. The average value of the F^ 2 tensor for this 
material was determined to be - 6.387 x 10 -4 (K.S.I.)" 2 . 

Even though the experimentally determined value of 
F 12 lies in the range - 5.933 x 10~ 4 to - 6.690 x 10" 4 (K.S.I.)” 2 
tension-compression tests were used instead of compression- 
compression tests, as is indicated in Figure 19, to obtain 
the best resolution of the F^ 2 tensor because the latter mode 
of testing was quite inconvenient to use with tubular speci- 
mens and this testing machine. Figures 23 and 24 show typical 
failure modes for the F-^ 2 test specimens. Other failure 
tests were conducted to verify this value of F 12 . Details 
of these results and analysis of the cubic int eraction terms 
can be found in Section 6 . 

F 16 ant ^ F 26 interaction Tensors 

The tests to determine the values of F. and F_ - 

16 26 

strength tensors have not yet been done. However, the optimal 
biaxial stress ratios for the four stress states and the 
attainable resolutions of the tensors have been calculated. 

These are plotted in Figures 25 to 28 and are summarized in 
Table III. 

However, it should be noted that since F g = 0, this 
would indicate that the sign of the shear stress does not 
affect the failure stress. Consequently, to remove this shear 
stress "sign" effect from the failure equation, all odd-order 
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terms in CT g should be set to zero. Thus it is proposed to 
let = F^g = 0 in our subsequent failure analyses. This 

will of course also apply to the . cubic terms involving odd- 
order ^Tg components. 

3.3 The Effect of Temperature on the Strength Tensor 

Components ' 

A preliminary test programme was undertaken to eval- 
uate the effects of varying ambient temperature from - 80°F 
to + 250°F on the orthotropic elastic constants and burst 
strength of 3-ply, 90° laminated tubes. Appendix B contains 
a detailed discussion of the test methods used together with 
the required analysis. A summary of the experimental re- 
sults is contained in Table IV. It should be noted that two 
test conditions were investigated using internal pressure 
loading; (1) free axial expansion of ends, and (2) zero axial 
displacement due to imposed external constraints. The actual 
burst pressures recorded in Table IV are for case (1) load- 
ing. 

In order to demonstrate how this test data can be used 
to estimate the variation in strength tensors with temperature, 
one must first consider the ability of the strength criterion 
to predict failure at room temperature where the strength com- 
ponents are already known. For a 0 = 90°, 3-ply laminated 
tube under internal pressure, 

<Tl = T y = pR/3t 
‘T 2 ~ O'x “ pR/6t 
^6 * 0 


( 13 ) 
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Neglecting the cubic terms 
can be written as. 


F 122' F 112 ' 


the strength, equation 


4 <r 2 2 ( F„ + f 22 F,J ) + z<r 2 ( F, = 1 

■4“ 2 

Substituting for the strength components (see Tables II and 
III), the predicted failure stress is < 3"" 2 = 3365 psi. This 
agrees quite well with the average measured value of 3442 psi 
(within 3%) shown in Table IV at T = 70°F. Thus, the assump- 
tion of F 122 and F 112 being of second order (for 6 = 90°) 

appears valid. The interesting feature of Eq. (14) is that 
if one neglects the F^, F and F 1 terras for the 0 = 90° 

configuration tests, then the resulting approximation yields 
a failure stress ( CT 2 ) 3246 psi. This value lies within 
4% of the predicted strength based on the inclusion of all 
terras... Consequently, one can now obtain an estimate of the 
temperature variation of the F 22 and F 2 terms from the follow- 
ing strength equation; 


*2 F 22 + F 2 <X 2 * 1 


(15) 


where we shall define the reference temperature (RT) strength 
tensor components as 


F 2rt 


1 

y 


1 

Y' 


F22 


RT 


Y ' - Y 

YY’ 

1_ 

YY 1 


(16) 


We shall further assume that at some arbitrary temperature (T) , 
these strength components can be estimated by 
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^2r)rn 

F 2 < T > - 7I§ , 


F_ 0 (T) = F22 RT 

*22 ~tt "." ■i , 2 


(17) 


(1+k) 

where Y(T) = (1 + k)Y 

Y 1 (T) = (1 + k) Y' 
k = a constant 

Using the strength equation (15) to predict failure at some 
temperature T, one has 

F 22 (T) T 2 2(T) + F 2 (T) <;r 2 (T) = 1 
or substituting for F 22 (T) and F 2 (T) from Eq. (17) gives, 


F 2 rt < ^/ T ) ~ F 2 2 Rr Cr 2(0= O 


(18) 


where = (1 + k) . 

Consequently, knowing the F 2rt and F 22rt values as well as 
the failure stresses (T) at various temperatures, one can 
solve for ot(T) . Thus one can construct plots of the varia- 
tion in the strength tensor components F 2 (T) and F 22 (T) as 
a function of temperature. 

Based on the data shown in Table IV for circular cyl 
indrical tubes with a mean radius R = 1.015" and a ply thick 
ness t = 0.010", the strength terms were then calculated for 
the various temperatures. These results are shown in Figure 
29. 


3.4 The Effect of Post Cure Time on the Strength Tensor 
Components 


Similar to the temperature study, an investigation 
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was made on the effects of varying post cure times (0 to 24 
hrs) on the measured orthotropic elastic constants and the 
burst pressure. A description of this work can be found in 
Appendix C, with a summary of the results listed in Table V. 
Since the test procedures and analysis are identical to the 
thermal study previously discussed, it is only necessary to 
consider the appropriate strength tensors F 2 and F 22 . 

Based on the post cure data, an average failure 

stress of = 3112 psi was obtained. This can be compared 
with the predicted (room temperature) value of 3365 psi for 
a 3 ply, 0 =90° configuration which agrees within 8%. 

Again, using the same argument posed in Section 3.3, the var- 
iation in and F 22 with post cure time was estimated and 
the results plotted in Fig. 30. 
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4. MATERIAL PROPERTY QIARACTERIZATIOiSr TESTS 


Although this report is primarily concerned with the evalu- 
ation of a strength criterion for the design of composite structures, 
it is also necessary to include the determination of the lamina 
orthotropic stiffness parameters. This arises because in the analysis 
of laminated configurations consisting of many plies oriented at 
arbitrary angles e relative to the structural axes, it is necessary 
to ocnvpute the individual ply stresses (aj, 02 , ci 6 or a , a , t ) 

x y ■ xy 

and strains. This will be discussed later in Section 7. Suffice 
it to say at this point that the plane stress distrubution in some 
lamina of orientation 0 (see Fig. 3) is related to the overall struc- 
tural strain matrix in the following way; 


where 


0 


X 


a 

= 

Y 


T 


L xyJ 

A *- 


On 

012 

OlG 

01 

ro 

1 — 1 

Q 22 

026 

Qei 

062 

066 

> cos 

20 + 


> cos 

20 + 

u 3 
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X 


e 


y 

j . 

-< 

• — 


(IS) 


0l2 = 021 “ U 4 ” U 3 COS 49 

Qi 6 = Q 61 =(U 2 /2)sin 29 + U 3 sin 40 (20) 

026 = 062 =(U 2 /2)sin 26 - U 3 sin 40 


066 = U 5 - U 3 cos 46 


and 


Ul = |(3 Qn + 3 Q 22 + 2 Q 12 + 4 Q 66 ) 
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°2 - " Q22) 

u 3 = |(Qn + Q22 - 2 Q 12 - 4 Q 66 ) (21) 

Ut, = £-(Qii + Q22 + 6 Q12 - 4 q 66 ) 

u 5 = ^(Qll + Q22 - 2 Q 12 + 4 Q 66 ) 

and 

Qll = E ll/(1 ~ vi 2 \> 2 i) 

Q 22 — E 22 /(l - vi 2 v 2 i) (22) 

Ql 2 ” Q2I “ V 21 E ll/d “ v 12 v 2l) 

Qie - Q 2 g = 0 for an orthotropic lamina 
Qs 6 = G i2 

Equations (19) to (22) are based on three major assumptions : 

(a) linear response; 

(b) the laminae exhibit identical structural strains under a given 
load system; 

(c) the individual laminae can be regarded as homogeneous orthotropic 
materials . 

Clearly the assumption of linear stress/strain behaviour requires 
experimental verification to define the limits for wliich this is 
true. Furthermore, one must also evaluate the actual orthotropic 
constants as defined by En, E 22 f v 12 (or v 2 i) and G 12 . 

The following sections describe the experimental methods 
used to determine the orthotropic constants and the results ob tain ed. 

4.1 EVaLUATICN OF OETUOTEOPIC ELASTIC CONSTANTS 

The material properties E n , E 22 , and G 12 were obtained 
from stress-strain data of tubes with 0 = 0° and 90° loaded in axial 
tension and tubes with 0 = 90° loaded in pure torsion respectively. 
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These material properties were calculated from basic relationships 
and are listed in Table VI. 

- Four other specimens were also characterized using the 

method outlined in Ref. 26. These specimens were all made up of 
four plies arid had the following fibre orientations: 0°, 0°, 0°, 0°; 

-30°, 30°, 30°, -30°; -60°, 60°, 60°, -60°; 90°, 90°, 90°, 90°. 

Three 350 ohm, 0.5 inch strain gauges were bonded onto each tube, 
oriented in the axial, circumferential and -45° directions. All of 
the strain gauges were positioned at the midlength of the specimens . 
Each tube was then potted in aluminum endplates using a rocm temper- 
ature curing epoxy. 

The characterization procedure consisted of applying 
separately to each tube pure axial canpression, internal pressure 
and pure torsion. By monitoring the load transducers and strain 
gauges on an X-Y plotter, a plot of load versus strain was obtained 
during each of the loading cycles for each of the strain gauges. 

Frcm these curves, the nine values of load/strain were calculated 
for each tube . These values along wdth the ply thickness and fibre 
orientation were the input data for the material properties charac- 
terization computer program listed in Ref. 26. This program cal- 
culated, among other things, the material properties E 22 , Gi 2 > 

\>2 1 and v ] 2 and these are also listed in Table VI along with the 
manufacturer ' s specifications . 

It can readily be seen by comparing the unidirectional 
tube data with the calculated values for the + 6 laminated tubes 
that for the E 11# E 2 2 and G^ 2 parameters, the results agree within 11%. 
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In the subsequent failure analyses, the unidirectional tube data 
will be used for these three orthotrqpic constants, while the major 
and minor Poisson ratio terms obtained from the + 0 tests will be 
employed . 

Of particular interest for stress calculations is the range 
of linear behaviour as shewn in Figs. 9, 14 and 15. From this data, 
the following limits should be considered on the use of Eq. (19) : 

- 88 ksi <_ ai _< 121 ksi 

- 6 ksi <_ c 2 1 . 3.3 ksi (23) 

1 06 ! 1. 3 ksi 

where 1 and 2 denote the fiber and transverse directions, respectively. 

4,2 THE EFFECT OF TEMPEfVmJEE ON THE ORTHOTBOPIC CONSTA^IS 

The effect of varying ambient temperature not only changes 
the strength parameters (as sho^vn previously in Section 3.3) but 
also the orthotropic stiffness coefficients. Appendix B contains 
a complete description of the tests performed over a temperature 
range of -80°F to +250°F. As far as the elastic constants are con- 
cerned, Figs. B.4 to B.10 illustrate the stress/strain curves obtained 
which are quite linear except near fracture in the high temperature 
range (T >_ 200°F) . The actual variations in the moduli Ej j and E 2 2 
are shown in Figs. B.ll to B.14 and summarized in Table IV. The 
major thermal effect is found in E 2 2 as expected since it reflects 
primarily the epoxy matrix response. Estimates of the temperature 
variation in the Poisson ratio terms were also made. Again, the 
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dominant change is in the \> 2 l parameter. 

Knowing both the change in strength and orthotropic 
stiffness parameters with temperature permits the designer to 
calculate the failure loads of composite structures at various 
thermal conditions. Although all parameters have not yet been 
determined as a function of temperature, the methodology outlined in 
this report should provide the means for obtaining the necessary 
data. 

4.3 THE EFFECT OF POST CUKE TIME ON THE OKTHOTBOPIC CONSTANTS 

As in the thermal investigation, the effect of varying 
the post cure time was studied to assess the change in the ortho- 
tropic constants. This aspect of the programme is discussed in 
Appendix C, with a summary of tire results contained in Table V. 

In general, it would appear that substantial variations in post cure 
time do not drastically alter either the stiffness or strength 
parameters. This is quite useful from a manufacturing viewpoint 
since over-curing of the material does not appear to be degrading. 
This of course would not hold true when one imposes too short a 
cure time. Other fabrication parameters such as cool-down rate 
should also be investigated to determine the tolerance limits in 
terms of their effect on mechanical properties. 

5. STRENGTH TESTS ON lAMENATED TUBES 


An extensive experimental investigation was undertaken to 
evaluate the strength of laminated 4 ply tubes having a symmetric 
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configuration (-0, +0, +0, -0) with respect to the mid-plane. Two 
different load conditions were employed , internal pressure and 
combined loading of internal pressure plus torsion. The major 
purpose of both series of tests was to obtain failure data as a 
function of fibre orientation to compare with predicted values based 
on the tensor polynomial strength criterion which is discussed in 
Section 6. 

5.1 BIAXIAL STRENGTH TESTS OF + 6 LAMINATED TUBES 


In the biaxial strength tests used to obtain the complex 
failure strengths, the tubes with the orientations listed below were 
loaded to failure by internal pressure. A total of fifteen four ply 
tubes were tested - three each of the following fibre orientations 


(see Table VII) ; 

0°, 0°, 0°, 0°; -30°, 30°, 30°, -30°; -45°, 45°, 45°, -45°; 

-60°, 60°, 60°, -60°; 90°, 90°, 90°, 90°. 


The same air operated hydraulic oil pump that was used in 
the F 12 tests was also used in this. set of tests (see Fig. 31) . The 
tube's internal pressure was again monitored by a pressure trans- 
ducer and its output was plotted on an X-Y recorder to obtain an 
accurate determination of the failure pressure. 

Fracture was evident by a "weeping" through the walls 
of the tube in which small droplets of oil appeared on the surface 
of the tube indicating that the matrix had failed. In the three 
tubes other than the circumferential and longitudinal wraps, this 
"weeping" was not localized at all, but rather, it occurred at many 
places around the circumference of the tube at the same time. The 
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internal pressure started to drop and could not be maintained even 
with an increased flew rate when the "weeping" started. Fracture 
was more pronounced in the 6=0° and e = 90° tubes as actual splits 
appeared in the matrix parallel to the fibre direction. A summary 
of the experimental results is contained in Table VIII. Later, in 
Section 6, a comparison between the predicted values and this data 
will be made. 

5.2 STRENGTH OF + 0 LAMUJATED TUBES SUBJECTED TO OTIEENAL PRESSURE AND 
TORSION 

Of the fifteen tubes manufactured for this investigation 
(see Table IX) , twelve ware tested to failure under combined torsion 
and internal pressure loading. Each tube was first subjected to 
a positive constant torque of 500 in. -lb. and then loaded to failure 
under internal pressure. Hie constant torque was applied using a 
Tinius Olsen universal testing machine modified for torsional loading. 
Internal pressure was generated using a specially designed hydraulic 
pump. In each case, the test specimen was completely filled with 
hydraulic fluid before being mounted into the testing machine, with 
care taken to bleed all trapped air frem the hydraulic lines. 

Hydraulic pressures were recorded using a Kistler Model 607A quartz 
pressure transducer with the output plotted on an X-Y recorder. 

The test specimens were in all cases judged to have failed 
when the first signs of hydraulic fluid could be seen on the outer 
surface of the tube. Accompanying these first signs of hydraulic 
fluid was an irrr, mediate drop in the recorded pressure. 

Figure 32 shows a photograph of the test apparatus used 
to accomplish the required combined loading. A summary of the failure 
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pressures can be found in Table X. Again, as with the biaxial 
strength tests, a comparison with predicted values will be done in 
Section 6 . 

6. EVALUATION OF THE TENSOR POLYNOMIAL FAILURE CRITERION 

The general form of the tensor polynomial failure criterion 
proposed by Wu is 


F. a . + F . a . o . 
ii iJ i 3 


a . 

i 


a . 


3 


°k 


+ 


( < 1 no failure 

= 1 failure 

> 1 exceeded failure 


(24) 


for i, j, k, = 1, 2, 3, ... 6. F . , F . . and F . .. are strength tensors 

i iJ 1D K 

of the 2nd, 4th and 6th rank, respectively. If one restricts the 
analysis to a plane stress state and considers only a cubic formul- 
ation as being a reasonable representation of the failure surface, 
then Eq. (24) reduces to 

F l°l + F 2 C 2 + F 6°G + F ll a l 2 + F 22°2 2 + F 66 a 6 2 

+ 2F i2°l°2 + 2F 16°l a 6 + 2F 26 a 2 cr 6 + ^ F 116 a l 2 °6 
+ 3F i26 cr l <J 2 cr 6 + 3F 112 (7 l 2cr 2 + ■ 3F 221 a l l7 2 2 + 3F 166 a l°6 2 
+ 3F226 a 2 2 °6 + 3F 2 6602°6 2 + F lll a l 3 + F 222 a 2 3 + F 666 a 6 3 = 1 (25) 



if it is further assumed that the material has sane form of symmetry 

(Ref. 6) such that F. . = F.. for i ^ j and F. = F.. , = F... - 

31 J 13k iki 31k 

Fjj,^ = = Fj_ Since it has also been shown (Ref. 6) that 

inclusion of the cubic terms F. . . (for i = 1. 2 and 6) is redundant. 

ni 

therefore they can be emitted. One other important simplification 
of Eq. (25) can also result if it can be experimentally determined 
that a lamina exhibits identical strength for both positive and 
negative shear. If this condition is satisfied, then all odd-order 
terms in can be set to zero to remove the shear stress "sign" 
dependence. Hence Eq. (25) reduces to 


F l°l + F 2 a 2 + F ll a l 2 + F 22 ci 2 Z + F 66 a 6 2 
+ 2F^ 201^2 "F -^ F 112^1 2 ®2 3 f 221*^2 2 ^1 

+ SFiggffiOg 2 + 3F 266 o 2 a 6 2 = 1 

Based on the experimental results (refer to Table I) used 
to evaluate the strength tensor .components (F^ and F__ ) described 
in Section 3, it was found (Fig. 15) -that the sign of the shear 
stress did not significantly affect the failure stress. Thus Eq. 

(26) is proposed as a failure criterion for the glass-epoxy material 
being studied. As shown earlier, the following strength tensor 
components have been determined experimentally; 
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Fj = -3.076 x 10" 3 (ksi) -1 


F 2 = 2.344 x KT 1 (ksi)” 1 


F 6 = 0 

F u = 9.398 x 10 -5 (ksi)" 2 (27) 

F 22 = 2.270 x lo- 2 (ksi) -2 
F 66 « 2.142 x 10 -2 (ksi) -2 
F 12 = -6.387 x l(T 4 (ksi) -2 

ESTIMATE OF CUBIC HJTEPACTION TERMS . (E^^) 

Since no direct experimental measurements were made with the 
intension of evaluating the cubic interaction terms, it is necessary 
at this stage to attempt to estimate their values based on the test 
data obtained. 

It can be shown that as 0 -> 60°, the shear stress term 
as becomes very small for internal pressure loading acting alone. 

However , since the failure pressure increases by over a factor of 
five compared to the 0° laminate (see Table VIII) , then the inter- 
action terms F 122 and F 112 must play a significant role for this 
60° orientation. On the other hand, since the quadratic form of the 
failure criterion provides accurate estimates of the strength for 



38 


the 0° and 90° configurations, then the cubic terms must have a small 
influence in these fiber angle ranges. Similarly, one can argue 
that for the internal pressure + torque tests, the major contributions 
frcrn the cubic terms should occur in the 40°-70° range. Consequently, 
to evaluate the four cubic terms F 112 , Fj 22 i F 166» F 2 6 6/ At was 
assumed that the other tensor components ■were known (i.e., see Eq. 

(27)) and the cubic polynomial (Eq. 26) was fitted to the test data 
at 45° and 60° for both the pressure and pressure-torsion cases. 

This yielded the following values for the cubic terms; 

F 112 = 1.6671 x 10 _lf (ksi)" 3 

Fi 22 = -1.0575 x 10“ 3 (ksi) -3 

f 

- " (28) 

F l66 = 3.9985 x 10- 4 (ksi)“ 3 

F 2 66 = -1.3726 * 10“ 3 (ksi)" 3 

One of the difficulties associated with a cubic formulation 
is that of ensuring a real positive root corresponding to a failure 
pressure. This of course depends on the value of the discriminant 
of the coefficients. For example, consider a general cubic polynomial 
of the form 

x 3 + bx 2 + cx + d = 0 (29) 


Hence, the discriminant can be defined by 
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A = 18 bed - 4 b 3 d + b 2 c 2 - 4c 3 - 27 d 2 
< 0 one real root, two canplex 
= 0 all roots real, two are equal 
> 0 three roots real and unequal 

Since the strength coefficients vary from the mean values defined by 
Eqs. (27) and (28) , it is quite possible for A to assume any of the 
three conditions stated above for a particular laminate configuration. 
This is readily apparent in Fig. 33 where the cubic failure poly- 
nomial has been plotted for the + 45° laminae. In one case (0 = -45°), 
two distinct positive roots were obtained differing by a factor of 
two. Hence, from a physical viewpoint, the lowest value should be 
taken. On the other hand, for 9 = +45°, the polynomial does not 
cross, but the "error” is very small. In other words, because of 
coefficient variations, one can consider the maximum as a "double' 
root", thus yielding a failure load. 

One other comparison that can be made is between the quad- 
ratic and cubic failure predictions. Figure 34 contains both poly- 
nomial solutions for the case of symmetric laminated tubes subjected 
to internal pressure. In this instance, failure of the + 6 and - 6 
laminae is assumed to occur at the same load. It is quite evident 
that the cubic terms are most significant in the fiber angle range 
of 40° to 70°. Although the quadratic solution provides the best 
estimates of the failure pressures at 0° and 90°, there is some 
discrepancy in the cubic solution near 90°. One other feature worth 
noting is the somewhat "clipped" response of the cubic prediction 
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near the peak. It might well be that small variations in the cubic 
parameters could result in a smoother distribution. In addition, 
no correction was made to the interaction parameter F 12 due to the 
presence of the F 12 2 and F 122 terms. Again, this nay lead to a 
different response. 

7. COMPARISON OF PREDICTED STRENGTH WITH TEST DATA 


This section of the report presents the analysis required 
to carpute both the principal lamina stresses and the failure stresses 
as a function of ply orientation. The purpose here is to test the 
validity of the cubic approximation of the tensor polynomial by 
comparing predicted failure loads with the test data contained in 
Section 5 . ~~ 

7.1' CALCUhATION OF IAMINA PRINCIPAL STRESSES 

In order to utilize the failure equation, one is first 
confronted with the calculation of the lamina principal stresses. 
Assuming linear elastic behaviour, the lamina inplane stresses 
corresponding to the structural axes (x,y) can be determined for any 
ply orientation (6) from Eqs. (19). This relation further assumes 
that for a given external load system, the structural strains define 
the laminae strains. Hence, based on known applied loads, the strain 
matrix can be determined frcm 


e 


7 

_1 

" 

N 

X 


A. . 


X 

E 




N 

Y 




y 

y 




N 





xy 





L J 



41 


where the load matrix represents forces/unit length and each of the 
terms is given by 


n 


A ij jj x Q ij (t k " t k-l ) 


where k - 1, 2 , ... n defines the laminae . For specially ortho- 
tropic construction (which exists for a symmetric laminate) , then 


( 32 ) 


Aie - A 26 = 0 


(33) 


Consequently, 


[A ij r ' 


A 22 

A-12 


A 

a 12 . 

An 

A 

. ~ 

0 

0 


1 

A 66 


where A = Ajj A 2 2 “ A 12 2 . Substituting Eqs. (31) and (34) into 
Eq. (19) yields the lamina stresses 


(34) 


O 

X 


Qii 

Ql2 

1 

Ql6 

0 

y 

rs 

Q 21 

Q22 

Q26 


0 

Qei 

L 

C>62 

Q66 


a 22 /a -a 12 /a 0 

-A 12 /A A u /A 0 

o 0 l/Aee 


N 


x 


N 


N 


xy 


(35) 


The final step in the analysis is the transformation from the struc- 
tural axes (x,y) to the principal lamina axes (1,2) for each ply 
having an orientation 0 with respect to the structural axes. This 
yields the in-plane principal lamina stresses, 



where m = cos 6 
n = sin 0 

7.2 CALCULATION OF FAILURE STRESSES FOR A LM-tENATED STRUCTURE 

Using the failure criterion defined by Eq. (26) where the 
sign of the shear stress (a 6 ) is assumed not to affect the strength 
prediction, simultaneous failure of all laminae can occur for certain 
load conditions (such as internal pressure for example) if the struc- 
ture is of special orthotropic construction (A 16 = A 2 6 = 0) . In 
general" however , for arbitrary loading (such as internal pressure 
and torsion for example) , the individual laminae will have differ- 
ent principal stress distributions even if the configuration is 
specially orthotropic. For these cases, each lamina must be 
analyzed separately and the failure criterion applied to determine 
which lamina (or pairs of laminae) will fail first. Subsequently, 
the stresses must be adjusted due to partial structural failure and 
the failure condition applied again. This multi-mode type of 
failure analysis should then provide an estimate of the overall 
failure loads. 

7.3 FAILURE PREDICTIONS FOR THE INTERNAL PRESSURE TESTS ON SYMMETRIC 
LAMINATED TUBES 

Using the mean values of the orthotropic elastic constants 


where 


m+e = 


m 2 n 2 
n 2 m 2 


+ 2 mi 


+ 2 inn 


+ mn + mn (m 2 - n 2 ) 


listed in Table Vi, 
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E n = 0.477 x 10 7 psi 
E 22 = 0.121 x 10 7 psi 
G 12 = 0.420 x 10 6 psi 
V 12 = 0.336 
v 2 i = 0.0855 

the orthotropic stiffness parameters can be calculated , 


(38) 


Gil ~ E n /(1 - v 12 v 2 i) “ 0.491 x lo 7 psi 

Q 22 = E 22 / (1 — vj 2 v 2 ^) ~ 0.125 x 10 7 psi 

Ql2 = Q 21 ~ V 21 E 2 1 / (1 - v 12 v 21 ) = 0.042 x 10 7 psi 

Ql6 = Q61 ~ 0 

Q 66 " G 12 = 0.042 X lO 7 psi 


(39) 


Hence the general orthotropic stiffness coefficients can be readily 
determined as a function of 6 using Eqs. (20) and (21). 

For the case of internal pressure loading with the ends 
free to expand axially, the generalized load matrix is given by 

1/2 

= pR 1 

0 

L 

Thus, substituting Eq. (40) and the appropriate Cl j , terms into 
Eg. (35) yields the lamina in-plane stresses which will take the 
following form; 





44 


a 


' 0. 1250* 

X 



a , 

m 

0.2500 

y 

t 


T I 


x 


e 

xy 



_ . 


where t = lamina thickness. It is assumed that the laminate is 
comprised of four equal thickness laminae. Hence the and stresses 
are identical for each ply orientation 0 and only the shear stress 
varies with 6. The final step in the analysis once t (6) is 
known is the multiplication by the transformation matrix [T] +g given 
by Eq. (37) . 

Based on this method of analysis, a computer programme was 
written to determine the lamina principal stresses and solve 
the cubic equation (26) defining the failure condition (see Appendix 
D) . The appropriate tensor strength components given by Eqs. (27) 
and (28) were used as the coefficients of Eq. (26) . The 
actual failure, pressures are recorded in Table VIII and a graphical 
comparison with the experimental values is presented in Fig. 35 as a 
function of ply orientation 0. Although the cubic terms were obtained 
from solutions of the polynomial matched to the 45° and 60° data, the 
remaining test results are also in reasonable agreement with the 
predicted response. Figure 36 presents the variation in principal 
failure stresses with fiber angle. These stresses can then be com- 
pared with the linear elastic limit values shown in Figs. 9, 14 and 
15. It can be seen that the shear stress (a 6 ) is only slightly in the 
nonlinear range and thus the assumption of linear elastic behaviour 
is acceptable. Moreover, it would appear that the failure mode 
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throughout the whole range of configurations was defined by a matrix 

tensile fracture. . ~ 

7.4 FAILURE PREDICTIONS FOR THE COMBINED LOADING TESTS (INTERNAL 
PRESSURE AND TORSION) ON SYMMETRIC LAf’-lINATED TUBES 

Although the tube configurations studied here are similar 
to those described previously (4 ply symmetric laminates) , the load 
condition has been complicated by the addition of a constant torque 
(- 500 in-lbs) to the internal pressure loads. For pure torsion (T) 
and constant shear stress across the laminate thickness t (= 4t ) , 
then the following equation holds , 


T = 2N X yA 


where A - enclosed area 


N xy ' Tt 




= 0.0765 kips/in 


(42) 


(43) 


Thus the generalized load matrix for this case is given by 


N 

X 


plV2 

N 

y 

= 

pR 

N xy 


.0765 


( 44 ) 



It can readily be shown that with the addition of the torque 
term, the principal stresses in the +0 and -9 laminae differ signifi- 
cantly. Consequently, a failure prediction must be> made for both sets 
of laminae at each fiber orientation to determine which laminae will 
fail first (see computer programme in Appendix d) . The results of 
these calculations are shown in Fig. 37 where both lower bound and 
median response curves have been indicated. Of seme interest is the 
fact that initial failure conditions can occur in different laminae, 
depending on the fiber orientation. The lower bound curve can be 
used as a conservative estimate of the overall failure load if one 
assumes that once a set of laminae have fractured, the remaining 
structure cannot withstand the resulting load and thus coincident 
failure results. However, one can expect that seme residual strength 
remains. even after initial matrix failure in the first set of laminae 
and coupled with the other laminae, the ultimate structural failure 
stress should be somewhat higher than the lower bound predictions. 
Consequently, the median response curve can be -used as a first approx- 
imation for comparison with the experimental failure loads. This 
was done in Fig. 38 where once again, general agreement was obtained. 

Based on this failure analysis, the predicted principal 
stresses are shown for both +9 laminae in Fig. 39. Again it can be 
seen that the assumption of linear elastic behaviour was reasonable 
and the failure mode was one of matrix tensile fracture. 

8. CONCmSIONS 

High quality laminated tubular specimens have been manu- 
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factored using a belt wrapper apparatus and subjected to the 
following load conditions: uniaxial tension and compression, pure 

torsion, internal pressure and combined torsion and internal pressure 
Failure loads were obtained to evaluate the principal tensor strength 
parameters (F^, F^; i = 1,2,6) and the quadratic interaction tem 
F 12 . It was then shown that a quadratic failure formulation was too 
conservative and a cubic polynomial can in fact be used to predict 
with reasonable confidence the failure loads for arbitrary fiber 
orientations. This was done on the basis of estimating four cubic 
parameters (F 12 2 / F i22* F i66/ F 266 ) using test data at 0 = 45° and 
60°. Consequently, it is highly desirable to not only obtain more 
test data for comparison purposes but refine the estimates of these 
cubic terms. In addition, various failure modes of the tubes were 
considered. The case of internal pressure loading of symmetric 
laminated tubes resulted in simultaneous failure of all laminae. 
However, the addition of torsion to the loading indicated that fail- 
ure of sane laminae should precede overall structural failure. 

Experiments were also performed to assess the effects of 
ambient temperature and post cure time on the orthotropic properties 
as well as on the failure parameters. Although only limited data 
were obtained, it was demonstrated ha/ these effects can be inter- 
preted and sane results were presented for two strength parameters. 
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TABLE I 


TENSILE, COMPRESSIVE, AMD SHEAR FRACTURE DATA 
SPECIMEN NO. ULTIMATE STRENGTH (P.S.I.) 

Longitudinal Tensile Strength Tests 


2a 

0* 

131,158 

2b 

0* 

111,004 

3a 

0° 

119,735 

3b 

0° 

121,336 

Mean 

Value 

120,303 


Longitudinal Compressive Strength Tests 


7a 

o°- 

33,784 

7b 

0 # 

89,491 

3a 

0* 

86,013 

9a 

0 ® 

86,885 

9b 

0° 

91,629 

13a 

0® 

88,676 

19a 

o’ 

87,091 

Mean 

Value 

83,081 


Transverse 

Tensile Strength Tests 


la 

90* 

ID 


2,980 

la 

90® 

(ii) 


3,520 

3a 

5a 

90® 

90® 

(D 

^°**'** uM 

3,204 

3,150 

5a 

90° 


3,333 

6b 

90° 

. — ■ 

3,418 

8b 

90* 



3,109 


3,246 


Mean Value 



TABLE I (cont*d«) 


TENSILE/ COMPRESSIVE AND SHEAR FRACTURE DATA 


SPECIMEN NO. 


ULTIMATE STRENGTH (P.S.I.) 


Transverse Compressive Strength Tests 


2a 

90° 

13,153 

3b 

90* 

14,035 

4b 

90* 

14,246 

7b 

90° 

12,213 

Sc 

90° 

13,319 

16b 

90* 

13,431 

17a 

90* 

14,114 

Mean 

Value 

13,574 


Shear 

2b 

Strength Tests 
90* 

6,737 

4a 

90*(I) 

6,777 

4a 

90*(II) 

6, $51 

5b 

90* 

6,S64 

6a 

90* 

6,765 

7a 

90* 

6,497 

17b 

90* 

7,336 

Mean 

Value 

6,332 




TABLE II 


STRENGTH TENSORS 


li 

-3.076 X 

10-3 

K.S.I.-l 

11 

CM 

2.344 x 

10-1 

K.S.I .” 1 

F6 = 

0,0 


K.S.I . -1 

F ll = 

9.39^ x 

10-5 

K.S.I .- 2 

f 22 = 

2.270 X 

10- 2 

K.S.I .- 2 

F 66 = 

2.142 X 

10-2 

K.S.I . “2 

|fi 2 |« 

1.460 X 

10-3 

K.S.I. -2 


1.419 X 

t 

O 

H 

K.S.I. “2 

M~ 

2.205 X 

lO - 2 

K.S.I . -2 



TABLE III 


TENSOR 
Fi2^-0. 5 
F 12^S °‘5 
F 16^0 
f 16^ 0 
f 26^ 0 
F 26^ 0 

SPECIMEN 

NUMBER 

10a 90* 

10b 90° 

11a 90° 

13b 90° 


INTERACTION STRENGTH TENSORS 


OPTIMAL STRESS STATE 

0*1 >0, 0*2*0 
0*i< 0, 0*2 ^ 0 

0'i>0, 0*5 <0 

<r 1 >o, cT6>o 
cr 2 ^o, <T6 >o 
<T 2 <£0, O*6<0 


RANGE OF OPTIMAL B 
-14.4^0i/<T 2 ^-15.6 
15.5^i/fi^ ^ 13.3 
-15 . 65^1/^ £--16.35 
16 . 35 * 1*1 /% > 15.65 
-O.S9>06/®2 >-1.0 
1 , 0? *£/02 > -0.S9 



F-12 TENSOR 

TESTS 




PREDICTED 

MEASURED 

ULTIMATE 

COMPUTED 

l/ 2 

l/ 2 

STRESS (K.S.I.) 

Fi p (K.S.I.-^ 

) 

-S.00 

-S.69 

Ti = 

10S.S 

1.S69 

x 10"4 



«*2 = 

-12. SS - 



-14.50 

-15.49 

0*1 = 

103.4 

-6.690 

x 10-4 ' 



• 

02 = 

-6.392 




-13 .SO 

-13.33 

<Ti = 

100.0 

-6.539 

X 

H 

O 

1 

>* 



cr 2 = 

-6.911 




-13.90 

-15.13 

0-1 = 

105.4 

-5.933 

x 10“4 j 



0*2 = -6.669 


Mean Value of F12 = -6.3S7 


x 10-4 . K.S.I. -2 



TABLE IV ELASTIC CONSTANTS AND TRANSVERSE BURST STRENGTH FOR SCOTCHPLY (1002) TUBES AT VARIOUS 

TEST TEMPERATURES 


TEST 

TUBE 

TEMPERA- 

DESIG- 

TURE 

NATION 

- 80 *F 

24b 90° 

- 75°F 

23c 90° 

- 45°F 

23b 90° 

- 20°F 

23a 90° 

- 20 °F 

22b 90° 

. - 20®F 

22c 90° 

70 °F 

12c 90* 

70°F 

13c 90° 

70 *F 

11c 90° 

158 a F 

19b 90° 

158°F 

16a 90° 

200®F 

19a 90° 

200°F 

! 14a 90° 

250°F 

14b 90® 

250 d F 

14c 90° 

! 250*F 

15a 90® 


E ll/0 

(psi) 


5 . 76x106 


5.86xl0 6 


6.01xl0 6 


6.40xl0 6 


6.35xl0 6 


6.0Sx1q6 


6.30xl0 6 


6.0Sxlo6 


e 22 / 3 

(psi) 


En 

(psi) 


e 22 

(psi) 


V12 




1,76x106 


1.39xl0 6 


l-.51xl0 6 


6,3 2x1q6 


5 . 04x1q6 


5.80x10 6 


1.73xl0 6 0.250 


1.22x10° I 0.721 


1.46xl0 6 0.356 


1.32x10° 


1.29xl0 6 


l.lSxlO 6 


1.15xl0 6 


o.98xio 6 


5 . 6 1x10 6 


5 • 66x10^ 


6.11xl0 6 


6.07X10 6 - 


5 . 7$x1q6 


1 . 26 xl 0 6 


1 '10..1 r\ 6 


1.22x10 


1.13x106 


I.IOxIO 6 


0.93xl0 6 


0.437 


0.520 


0.495 


0.495 


0.552 


0 . 5 7xl0 e 


5.S7xl0 6 


0.66x1q6 I 5.91xl0 6 


0.53xl0 6 0.875 


0. 65 x 10 6 0.511 



6.23x106 1.59x1q6 5.86x10 6 I 1.50xl0 6 | 0.482 0.123 


0.068 


0.173 


0.089 


0.098 


0.112 


0.092 


0.090 


0.089 


0.079 


O.O 56 


OVburst 

A* 

(psi) 

(Mode 1 ) 







4399 


4186 


4110 


4200 


4224 


3411 


3369 


3545 


3216 


2876 


2300 


2450 


1405 


1280 


1325 












































































































TABLE V ELASTIC CONSTANTS AND TRANSVERSE BURST STRENGTH FOR SCOTCHPLT (1002) TUBES HAVING 

DIFFERENT LENGTHS OF POST CURE 


LENGTH OF 
POST CURE 

TUBE 

DESIG- 

NATION 

Ell/3 

(psi) 

E22/3 

(psi) 



Ida 90 ° 

5 . 896 x 106 

1.307xlo6 



18 b 90 * 

6.049x106 

1.417xl0 6 



18 c 90 * 





5.670x106 


7 

7 


2 


HRS 

HRS 




21 a 90 * 


21b 90° 

6.054xl0 6 

1.329xl0 6 

5.599x106 

2lc 90* 

5.855x106 

1.260x106 

5.543x106 

11c 90° 

6.013xl0 6 

1.297xl0 6 

5.663x10 6 


12 c 90 ° 


13c 90° 

5.865x106 

1.319x106 

5.612x106 

20a 90* 

5.978x106 

1.283xl0 6 

5o743xl0 6 

20b 90° 

5.908x106 

1»273x10 6 

5 O 426xlo6 


E 22 

(psi) 

1.239xl0 6 

1.328xl0 6 

1.229xl0 6 

1 . 193 xl 0 6 

1.222x10^. 

1.262xl0 6 

1.233x10^ 

1 . 169 x 10 ^ 





0.4374 

0.4280 


0,6152 


V 21 

^2 

(psi) 

(Mode l) 

0.1071 

3058 

0.1212 

2948 


3329 


3523 

0.1285 

3490 

0.1071 

3279 

0.1122 

3545 


3411 

0.0984 

3369 

0.0982 

3653 

0.1326 

3972 













































































TABLE YX 


MATERIAL PROPERTIES 


SPECIMEN 

^11 


E22 

G12 

V 12 

^21 

NUMBER 

(106 P.S. 

I.) 

(106 P.S. I.) 

(106 P.S. I. 

.) _ 


Calculated 

Directly From 

Stress-Strain 

Curves 



2b 0° 

4.74 






0 

O 

4.79 






Mean Value 

4.765 






3a 90° 



1.245 




5a 90® (I) 



1.157 




5a 90*(II) 



1.234 




6b 90° 



1.157 




Mean Value 



1.211 




2b 90® 




0.431 



4a 90*(I) 




0.403 



7a 90® 




0.420 



Mean Value 




0.420 



Calculated 

Using the 

Characterization 

X 

Computer Program 


lc 0° 

4.945 


1.170 

0.4333 

0.3043 

0.0720 

6b 0* 

4.676 


1.226 

0*4645 • 

0.2327 

0.07412 

2b + 30° 

5.107 


1.426 

0.4574 

0.2676 

0.07473 

7b + 30® 

5.447 


1.174 

O .5216 

0.4165 

0.03973 

4b + 60 ® 

5.031 


1.645 

0.4032 

0.2230 

0.07221 

9b + 60® 

5.125 


1.420 

0.4733 

0.2377 

0.0659 

5c 90° 

4.669 


1.243 

0.4215 

O .4074 

0.1039 

10b 90° 

4.376 


1.416 

0 . 4344 

0.4313 

0.1252 

Mean Values 

1 4.991 


1.341 

0.465 

0.3213 

0.03536 



TABLE VI (cont'd) 


MATERIAL PROPERTIES 


SPECIMEN Ell E22, Gl2 V 12 V21 

NUMBER (1QQ P.S.I.) (IQQ P.3.1.) (IQQ P.S.I.) 


Manufacturer's Values 


5.7 


1.4 



TABLE VII 


GEOMETRY 

OF TUBES USED IN BIAXIAL 

STRENGTH TESTS 

TUBE 

R 

7 

DESIGNATION 

(IN) 

(IN) 

6a 0 

1.017 

.0400 

6b 0 

1.017 

.0406 

6c 0 

1.017 

i 

.0407 

7a + 30° 

1.018 

.0397 

7b 7 30° 

1.018 

.0408 

7c + 30 

1.018 

. 0405 

8a + 45 

1.021 

.0406 

8b + 45 

1.021 

.0408 

8c + 45 

1.020 

.0397 

9a + 60° 

1.021 

.0401 

9b + 60° 

1.021 

.0410 

9c + 60° 

1.021 

. 0404 

10a 90° 

1.021 

.0410 

10b 90° 

1.021 

.0400 

10c 90° 

1.021 

.0404 



TABLE VIII 


BIAXIAL 

STRENGTH TEST DATA* 

SPECIMEN 

MEASURED - BURST 

NUMBER 

PRESSURE (P.S.I 

6a 0° 

131.0 

6b 0° 

135.0 

6c 0° 

108.0 

Mean Value 

124.7 

7a + 30° 

241.0 

7b + 30° 

260,0 

7c + 30° 

240.0 

Mean Value 

247.0 

'8a + 45° 

557.5 

8b + 45° 

582.0 

8c + 45° 

630.0 

Mean Value 

589.8 

9a + 60° 

675.0 

9b + 60° 

643.0 

9c +60° 

665.0 

Mean Value 

661.0 

10a 90° 

250.0 

10b 90° 

242.0 

10c 90° 

289.0 

Mean Value 

260.3 


*4 Ply Cylinders C-6, + 6, +6, -9) 



TABLE IX 


GEOMETRY OF LAMINATED TUBES FOR COMBINED LOADING TESTS 


TUBE 

DESIGNATION 

1 

0 

2 

3 

4 * 

L 

in. 

R 

in. 

t 

in . 

la 0° 

0° 

0° 

0° 

0° 

9.5 

1.016 

0.0399 

lb 0° 

0° 

0° 

0° 

0° 

9.5 

1.016 

0.0404 

1c 0° 

0° 

0° 

0° 

0° 

9.5 

1.016 

0.0401 

2a + 30° 

-30° 

30° 

30° 

-30° 

9.5 

1.019 

0.0408 

2b + 30° 

-30° 

30° 

30° 

-30° 

9.5 

1.020 

0.0410 

2c + 30° 

-30° 

30° 

30° 

-30° 

9.5 

1.019 

0.0406 

3a .+ 45° 

0 

<3- 

1 

45° 

45° 

-45° 

9.5 

1.022 

0.0409 

3b + 45° 

-45° 

45° 

45° 

-45° 

9.5 

1.022 

0.0410 

3c + 45° 

-45° 

45° 

45° 

-45° 

9.5 

1.022 

0.0404 

4a + 60° 

-60° 

60° 

60° 

-60° 

9.5 

1.022 

0.0405 

4b + 60° 

-60° 

60° 

60° 

-60° 

9.5 

1.022 

0.0414 

4c + 60° 

-60 

60° 

60° 

O 

o 

*£> 

1 

9.5 

1.022 

0.0410 

5a 90° 

90° 

o 

0 

90° 

90° 

9.5 

1.022 

0.0410 

5b 90° 

0 

o 

O'* 

90° 

90° 

90° 

9.5 

1.022 

0.0417 

5c 90° 

90° 

90° 

90° 

‘90° 

9.5 

1.021 

0.0410 


*lst lamina is innermost lamina 




TABLE X 


FAILURE PRESSURES UNDER COMBINED LOADING 
OF CONSTANT TORQUE AND INCREASING INTERNAL PRESSURE* 



FAILURE 

TUBE 

PRESSURE 

DESIGNATION 

(p . s, 1. ) 

la 0° 

112 

lb 0° 

145 

2b + 30° 

245 

2c + 30° 

276 

3a + 45° 

623 

3b + 45° 

597 

4a + 60° 

675 

4c + 60° 

597 

5a 90° 

276 

5b 90° 

256 

5c 90° 

.273 


*4 ply cylinders (-, + , -h, -) 




FIG. I, BELT WRAPPING APPARATUS 









ot 


FIG. 2 . SPECIMEN FABRICATION 




FIG. 4 SPECIMEN IN VACUUM BAG 
READY FOR CURING 



FIG. 5 APPARATUS FOR CUTTING SPECIMENS 




FIG. 6 LONGITUDINAL TENSILE 
STRENGTH TEST SET-UP 


* — 





FIG. 7 ONSET OF FAILURE 



FIG. 8 



COMPLETE FAILURE 


AXIAL STRESS (RSI) 



FIG.. 9 STRESS - STRAIN CURVES FOR SCOTCHPLY (1002) 
LOADED PARALLEL TO FIBER DIRECTION 




FIG. 10 LONGITUDINAL COMPRESSIVE STRENGTH TEST SPECIMEN 



FIG. II AFTER FAILURE 



FIG. 12 TRANSVERSE TENSILE FAILURE 



FIG. 13 TRANSVERSE COMPRESSIVE FAILURE 



FIG. 1 4 STRESS -STRAIN CURVES FOR SCOTCHPLY (1002) 
LOADED PERPENDICULAR TO FIBER DIRECTION 





FIG. 45 STRESS-STRAIN CURVES FOR SCOTCHPLY (1002) 
IN PURE SHEAR 






FIG. 16 SHEAR STRENGTH TEST SET-UP 



FIG. 17 SHEAR FAILURE OF A SPECIMEN 
LOADED UNDER PURE TORSION 
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FIG. Zl TEST SET-UP FOR THE DETERMINATION OF 
THE F , 2 INTERACTION TENSOR 


COMPRESSIVE LOAD (X I0' 3 VOLTS) 






r 
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FIG. 23 SPECIMEN FAILURE IN F, 2 TESTING 







FIG. 24 FAILED SPECIMENS FROM THE F (2 
TESTING 



FIG OPTIMUM RIAXIAI RATIO R, . FOR DETERMIMATION OF F,* 



FiG. 26 ATTAINABLE RESOLUTION OF F (6 FOR OPTIMAL RATIOS, B = -§!- 

€ 




FIG. 27 OPTIMUM BIAXIAL RATIO B FOR DETERMINATION OF Fas 

2o 
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FIG. 28 ATTAINABLE RESOLUTION OF F 26 FOR OPTIMAL RATIOS, B*- 2 ^ 

26 &Z 






FIG. 31 INTERNAL PRESSURE TEST SET-UP 



~ 

FIG. 32 COMBINED LOADING OF TORSION AND INTERNAL PRESSURE 


original TAGEIS 
OF POOR QUALITY! 


FIG. 33 PLOT OF CUBIC FAILURE POLYNOMIAL TO DETERMINE ROOTS 



FIG. 34 COMPARISON OF QUADRATIC AND CUBIC TENSOR POLYNOMIAL 
STRENGTH PREDICTIONS (ALL ODD-ORDER TERMS IN 
CTq SET TO ZERO) FOR SYMMETRIC LAMINATED TUBES 
UNDER INTERNAL PRESSURE 





101 (DEG.) 


FIG. 36 VARIATION IN PRINCIPAL FAILURE STRESSES WITH FIBER 
ORIENTATION FOR SYMMETRIC LAMINATED TUBES UNDER 
INTERNAL PRESSURE 



-6 




FIG. 37 COMPARISON OF PREDICTED FAILURE LOADS FOR ±9 LAMINATED 
TUBES UNDER INTERNAL PRESSURE WITH TORQUE USING CUBIC 
FAILURE TENSOR 




FIG. 38 COMPARISON OF MEAN CUBIC TENSOR POLYNOMIAL FAILURE 
PREDICTION WITH EXPERIMENT FOR SYMMETRIC LAMINATED 
TUBES UNDER INTERNAL PRESSURE AND CONSTANT TORQUE 



1 0| (DEG.) 


FIG. 39 PLOT OF THE PRINCIPAL MEMBRANE STRESSES AT 
FAILURE FOR SYMMETRIC LAMINATED TUBES UNDER 
INTERNAL PRESSURE AND CONSTANT TORQUE 



APPENDIX A 


DETERMINATION OF STRESSES IN FIBRE-REINFORCED CYLINDERS 
UNDER COMBINED LOADING 


A,1 Anisotropic Shell Theory 
1« Equilibrium Equations 

Let the mid-surface of the circular cylindrical shell be 
the reference surface with x, 6, and z measured with respect to that 
surface in the axial, circumferential, and radial directions respec- 
tively (Fig # A,l). The components of displacement u, v, and w of 
a point on the mid-surface are the displacements in the x, 0, and 
z directions due to the externally applied loadings. The thickness, 
h, of the tube is assumed to be uniform. Differentiation is denoted 
by a comma. 

For/ a cylindrical element, the strains are given in the 
following form: 







€ x - + Z K X 

— 6© + z[l — z/R] K 


(Al) 


^x© - Cl + z 2 /2R 2 ] €x® + z[l ~ z/2R] K x ^ 


where 


€x° = u >x 

= (1/R) (v,*+ w) 

6x« = (u, d ./R) + v, x 
K x = -w >xx 

K 0 = -(1/R 2 ) (w ># 0. + w) 

K x = -(2/R) [w, xd + (u,*/2R) - (v,x/2)] 


(A2) 


are the mid-plane strains and curvatures. 


Under axisymmetric loading differentiation with respect 
to the e coordinate is assumed to be zero. Substituting this 

C 

and equation (A2) into equation (Al) gives the strain-dis- 
placement relations in the following form: 


&x “ u >x ~ z w »xx 

€ © = ^jl_2 + z 2 

R R R 2 


Y® 


,= V 




1 + z 



(A3) 


The force and moment resultants are related to the 
middle-plane strains and curvatures in the following manner 



(Ref. 25): 


N x 


A 11 

a 12 

D i6 /2R 2 • 

D n/R 

0 


N © 


Al2 

A22 

D 2 6 / 2 R 2 

0 

-D22/R ■ 

-D26/2R 

Nxe 


0 

0 

a 66 +d 66/ 2r2 

D16/R 

0 

D66/2R 

M x 


Dn/R 

Di 2 /R 

D16/R 

Dll 

d 12 

Dl6 

M* 


0 

0 

0 

Dl2 

d 22 

D26 



D 16 /R 

D26A 

D66/R 

Dl6 

D26 

D66 

where 

(Aij, Djj,) -J" 

h/2 

Q-m (l,z 2 

-h/2 

)dz 

(i,j = 1 

,2,6) 



The equilibrium equations are obtained by considering 
the forces and moments acting on an infinitesimal element of 
the mid-plane of the shell (Fig.A>2). These are 


N x ,x = 0 


N xo>x + 1 M x«,x = 0 
R 


(A6) 


M x,xx. - Ne + P = 0 
R 


where p is the internal pressure applied to the tube. 

2, The Governing Differential Equation 

The first part of equation (A 6) gives, upon integration, 

that 

N x = constant = N 0 (A 7) 

by substituting the constitutive equations (A ; +) into the 
second and third parts of equation ^6) and using equation 



(A7) the following differential equation in w can be derived. 


Vl_w>xxxx + X 2 w >xx + X3W - P " A 12 A 66 No + Ai2^l6"T* (A^) 

R 2 R 2 RD 2 TTR 2 D 

where = ("dh - 1 ^D-^ 2 A55 - 4A^5 D16 Du + 4 A H Dl6 2 ) 

[ * 

X 2 = 2 (a66 A 12 Dll “ 2A 12 d 16 a 16^ % 

D 

X3 = A 22 - Ai 2 2 A66 

D 

A 22 = A 22 + P 22 
R 2 

a 66 = A 66 + 31*66 

R 2 

a 16 = »16 • 

R 2 

D = An A66 - ’ A^ 2 

3 . Solution of the Governing Differential Equation 

The differential equation (A&) can be solved by rewriting 
it in the form: 


w,xxxx 


+ X 2 W,XX + b w 

r 2 Xi r 2 Xi 



- A12 A 66 No + A 12 A16 




(A9) 


RD 


2 TTR 2 D / 


The solution to this equation is composed of a particular 



solution, w c , and the solution to the homogeneous equation, 
wji* For this case, the particular solution is 

w 0 = R 2 / P - A12 A66 N 0 + A12 A16 T 
5 ^ \ RD 2 TfR 2 D 

and the solution of the homogeneous differential equation is 
of the form 

w^ = k e mx (All ) 

where k and m are either real or complex constants. Sub- 
stituting equation (All) into the homogeneous differential 
equation and solving for m gives 

m 2 = -b +\/b 2 -4d (Al2) 

2 

where b = X2 and d = X3 
R 2 Xi R 2 Xl 

There are three possible cases that could be solved for, 
namely, 4d:>b^, 4d = b , and 4d<rb , however, only the first 
case will be solved in detail in- this paper. 

When 4d>b 2 , than equation (A 12 ) becomes 

m 2 = -|(-b + l\Jkd-b^) (A 13 ) 

= S + 1 0 




(A 10 ) 


where i =\^1 



Therefore, m can be written as 



m = ot _+ i 

(414) 

and 

m2 = cs.2 _g2 + 

(A15) 

so that 

© = «**2 _ /3 2 = _ b / 2 

(A16) 

and 

= 2e{/S= thA-d-b 2 

(AL7) 

Solving 

equation (Al6) and (A17) gives that 



* = i (2i/d - b)^ 

(A13) 

and 

= •§• ( 2/d* + b )^ 

(A19) 


Substituting these into equation (A14) and then into equation 
(AH) gives the following expression for the solution of the 
homogeneous differential equation 


^ (a cos jS x + B sin^x) +e _(4:x: (C cos/Sx + D sin/3x) (A20) 


and w = Wh + w Q (A21) 

The constants A, B, C, and D can be solved for through the 
use of the boundary conditions. 

Once the expression for w is known, its derivations can 
be calculated quite easily. In order to determine the stresses 
and strains in the tube, the only other quantities needed are 
u f x an d v yX , which can be calculated from the following 
relations 


u,x = 


^66 po ~ A l6 T 

[ A66 21T 


R 


- A 12 w + w jXX /Dn - SA-j^D^X 

R \ R RA66 J 


(A 22) 


D 


R 



04 23 ) 


V,x = A 11 “ A l 6 No + a 12 a i5 w + w >xx I2D15D-A15D11A66+2A16 2 D l 6 

2 TTRD D RD I RDA^b 

The stresses in each ply are calculated from Hooke’s law 


for plane stress, 


f— 

<Tx 

k 

Q 11 

Ql2 

§i6 

k 

~ -1 
^X 

k 




Ql2 

^22 

Q 26 



(A 24) 



Ql6 

Q 26 

Q66 


€. xd 



where Qij is the lamina stiffness matrix referred to an arbi- 
trary set of axes. 


A. 2 Conclusions 

In this paper IS presented a method derived from the 
extension of Flugge’s shell theory to anisotropic materials 
by which the stress variations along the length as well as 
across the thickness of a symmetrically laminated orthotropic 
cylinder can be calculated, A computer program has been 
written to perform this task and is listed in the following 
section. The results for a four ply circumferentially 
wrapped cylinder loaded under a combination of axial com- . 

i 

pression and internal pressure are also listed and plotted 
there. As can be readily seen from these results, the effects 
of the end contraints are confined to relatively small regions 
near the ends of the cylinder. 

This analysis allows the researcher to more optimally 
design composite cylinders so that the test section is suf- 
ficiently far from the ends of~the cylinder where the large 



bending stresses are induced- and the thickness of the tube 
can be chosen so that the stress variation across it remains 
very small. Also, other types of end attachments (boundary 
conditions) could be designed that might relieve these large 
bending stresses caused by rigidly clamping the ends of the 
cylinder. 


A .-3 Computer Program to Calculate Shell Stresses 


Purpose: To calculate (Tx, CFV , cr x ^ as functions of x, the 

position along the length of the tube, and z, the 
position across the thickness, for circumferentially 
wrapped tubes. 

Method: Calculate the mid-plane strains and curvatures and 

then determine the stresses from the Hooke 1 s law 
relations. 

Input: The data card contains the loading applied to the 

tube. 


Output : 



N 0 

1 

P 

~T~ 

T 

“1 

0 


10 

■ 

20 


30 


Position along the length of the tube, x (tube end 
at x=0); position across the thickness of the tube, 
z (mid-plane reference surface at z=0); applied 
loadings, N 0 , P»T; stresses Gr x ®; radial 
displacement and derivatives, w, w jX » w,xx« 


PAGE 1 ' ' ' ‘ : 

// JOB 0001 

LOG DRIVE CART SPEC CBRT AVAIL PHY DRIVE 
OOOO 0001 0001 0000 

V2 Mil ACTUAL 8K CONFIG 8K 


// FOR 

♦LIST SOURCE PROGRAM 
♦ONE WORD INTEGERS 
♦I OCS(1132PRINTER. PLOTTER* CARD ) 

C 

C THIS PROGRAM CALCULATES AS 

LENGTH OF THE TUBE) AND AS 
THICKNESS OF THE TUBE) THE 
THE SHEARING STRESS RESULTING 
FORCE* INTERNAL PRESSURE* AND 


C 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


A 

A 


FUNCTION OF X ( THE POSITION 
FUNCTION OF 2 ( THE POSITION 
AXIAL STRESS* THE RADIAL STRESS* 
FROM ANY GIVEN COMB INAT I OB OF 
TORQUE 


ALONG THE 
ACROSS THE 
AND 
AXIAL 


I 


REAL N0.H.L1 *L2»L3*LN*K< 4 ) 

DIMENSION Q(3»3)*E(3).STRl3) *COF (4*4 ) *W0M(4) *N1 (4) »N2 (4) 
PI *3.14159 . W 

R = 1 *020 
H=0o040 
Z=H/2 • 0 


$ri 4 LI L " Sri E LL 




THE NEXT FIVE CARDS GIVE THE ELASTIC CONSTANTS OF THE COMPOSITE 


E11=0.4765E7 
E22=0.1211E7 
Gl2 = 0 »4230E6 
V12= 0*3363 
V21=O.085S 


VV = 1 • 0-V 12*V21 

THE Q-MATRIX IS THE STIFFNESS MATRIX 


G( 1 * 1 ) =E22 /VV 
0(2*2 )=E11/VV 
Q { 1 * 2 ) =E 1 1 #V2 1 /VV 
0(2*1 ) =0 ( 1*2 J 
0(3*3 )=G12 
0(1.3 )*0.0 
Q(2.3 )=0*0 
0(3*1) =0(1 *3) 

Q( 3*2 ) =Q ( 2 *3 ) 
A11=0 ( 1 ♦ 1 >*H 
A1 2=0 ( 1 *2 )*H 
A 16=0 ( 1 » 3 ) *H 
A22=Q ( 2 ♦ 2 ) 
A26=Q(2*3)*H 
A66=0 ( 3*3 )#H 


JV.I Cl'! i (? i EtL PRO rir. I 5 TIE S 


■ 


-vt 


D11 = Q (1*1) *0/12*0 
D12=Q ( 1 *2 )*H**3» 0/12*0 
D2 2=0 ( 2 *2 ) ♦H**3* 0/12*0 
D16=Q (1*3 >*H**3«0/12«0 
D26=Q(2»3 )*H**3 *0/12*0 
D66=Q(3*3 )*H*«3. 0/12*0 






rtlPTCt K 


cut. me 


A.t 


rr>fiTR ix 
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TH l 
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THE NEXT 11 CARDS GIVE THE CONSTANTS USED IN THE D! FF » EON® 

F22=A22+D22/R**2«0 

F66=A66+3»0*D66/R#«2.0 

F16=D16/R**2.0 " 

D=A11*F66-F16*#2.0 

L1“D1 1“(DU**2.0*F66-4.0*011*016*F16+4.0*D16*#2.0*AU ) /<D*R**2«0) 
L2=2oO#(F66*A12*Dll-2.0*A12*F16*D16> /D 
L3=F22-A12**2.0*F66/D 
Cl=SQRT<L3/IR**2®0*Ll) ) 

C2=L2/(R**2»0*L1 ) 

D2=4.0#C1«*2.0 

B2=C2**2.0 

WRITE <3*500)B2»D2 

500 FORMAT I 1H »'B2= • *E10.4. 5X ♦ >402= ’»E10.4»///) 
ALF=5GRT(2oO*Cl-C2)/2«Q 

9ET=SQRT <2«0*Cl+C2)/2«0 

LN = 5 « 0 

C0=CC5<BET*LN) 

S I =S I N ( BET*LN ) 

EX=EXP<ALF*LN) 

EXM=EXP(-ALF#LN) 

THE COF-MATRIX IS FOUND BY APPLYING THE 6.Co‘S TO THE SOLUTION OF 
THE DIFF. EON. FOR W. IT IS THE MATRIX USED TO DETERMINE THE 
CONSTANTS A*B*C*D* 

THESE ARE FOR THE RIGIDLY CLAMPED CASE 

COF (1*11=1*0 
COF ( 1 »2 ) =0.0 
COF U *3) = 1.0 
COF ( 1 *4 ) =0.0 
COF ( 2 ♦ 1 ) =ALF 

COF (2*2) =BET ~ 

COF ( 2 » 3 ) =-ALF 
COF I 2 »4 ) =BET 
COF 1 3 » 1 ) =CO#EX 
COF ( 3 » 2 ) =S I #EX 
COF { 3 » 3 I =CO#£XM 
COF ( 3 *4 ) =S I *EXM " 

C0F<4«1 >=(ALF»CO-BET*SI) *EX 
COF (4 .2 )= (BET*C0+ALF*5I > *EX 
C0F(4»3)=(-ALF*C0-BET*SI > *EXM 
C0F(4»4)=(BET *CO-ALF*SI ) *EXM 
WRITE(3»600>C0F 

600 FORMAT ( 4E14.4 » / ) ' 


NOW FIND THE INVERSE OF COF TO FIND A»8*C*D 


CALL MINV(C0F»4.DET*N1»N2 ) 

WRITE (3»700)COF 
- 700 FORMAT UE14.4*/ ) 

READ! 2*100) NO»P*T 

100 FORMAT ( 3F10. 2 ) - 

CALL SCALF(6e0»0*10*0oG»-50.0) 
CALL FGR ID< 1 »0o 0 »-50a04 1 « 0* 100) 
CALL FGRIDtO«0*0»0.0«OolO»25 ) 
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CALL FCHAP(2*55.-1«0«0#2090.20.0.0) 

WRITE ! 7s 10) 

10 FORMAT ( 'X ( IN. ) ' ) 

CALL FCHAR(-0. 15.10*0*0. 20.0. 20.Pl/2.0) 

WRITE!7»11) 

11 FORMAT!' STRESS X (K«S«I< )') 

AX=0 . 50 

DO 12 1 = 1.5 

CALL FCHAR! AX-0 *05/6.0 .-2.50.0.10 .0*10. 3.0*P 1/2.0) 

WRITE ( 7 9 13 ) AX 
13 FORMAT (F3*l> 

AX=AX+0*50 

12 CONTINUE 
AV=-50.0 

DO 1 A 1=1.11 

CALL FCHAR (-0. 7/6*0. AY-0. 50.0*10.0.10 .0*0) 

WRITE ( 7»15)AY 
15 FORMAT ( F6 .1 ) 

AY=AY+10»0 
1 A CONTINUE 

CALL FCHAR (1.2 5 *-45*0.0. 20.0.20.0.0) 

WRITE! 7.81 ) 

81 FORMAT ( 'STRESS X AS A FUNCTION OF X FOR') 

CALL FCHAR! 1.25.-48*0.0. 20.0.20.0.0 ) 

WRITE < 7 .82 ) 

82 FORMAT! 'RIGIDLY CLAMPED ENDS') 

CALL FCHAR! 1.00.-51.0.0. 20.0.20.0*0) 

WRITE (7.83)N0.P»T 

83 FORMAT! »N0=* »F8*2. » LB./IN*.'.' P='.F7*1.' P.S.I..**' T='.F4*1.', ♦' I 

IN. -LB* ' ) ' ' ■ 

WO=R**2«0atP-1.0»< A12*F66*N0) /<R*D)+! A12*F16*T ) / l 2 *0*P I ) )-/L3/L3 
WOM(l)=-WO i;';. ' 

WOM ( 2 ) = 0*0 ' ; 

WOM ( 3 ) =-W0 * ' 

WOM (4) =0.0 r 

DO 3 11=1.4. f 

K ( 1 1 ) =0.0 

DO 3 J J® 1 ♦ 4 I \ 

K 1 II ) =K ( 1 1 ) +COF ! I I . JJ) #WOM ( J J ) — — 

3 CONTINUE J •• ' j 

WR ! TE ( 3 » 800 ) K f 

BOO FORMAT !4E14. 4) 

DO 21 L= 1 *3 

WR I TE ! 3 » 200 ) - 

200 FORMAT !1H1.3X»'X' *5X.'Z' . 5X . * NO ' »7X » * P * *9X » 1 T ' .5X s' STRESS X* .SX.'C.'s 
7TRESS Y' »5X . 'STRESS XY ' * 8 X > » W * ♦ 1 2X . ' WX 1 . 1 IX » ' WXX * . / ) 

X=0.0 ! 

DO 1 N= 1 . 51 }. 

CS*COS!BET*X> . ! 

SN=SI N ( 8ET*X ) } 

W = < K ! 1 )*CS+K!2>*SN)*EXP< ALF*X)+!K<3 )*CS+K!4)*SN >*EXP! ~ALF#X)+WO ) 
WX=< (K! 1 )*ALF+K! 2 ) *BET ) *CS+ ( K < 2 ) *ALF-K < 1 )*BET )#SN>*EXP (ALF*XH-< !X!(K( 
84)*BET-K(3)*ALF)*CS~(K!3 ) *BET+K ( 4 ) *ALF ) *SN )#EXP I -ALF*X > [ 

WXX=< <K! 1 )*ALF»*2.0 + 2.0*K; <2 >*ALF#EET-J<{ l)*6£T#*2.0)*CS-MK!2>ftALF*.F*# 
12 • 0-2 . 0*K ( 1 )'*ALF*BET-K< 2 ) *BET**2 .0 ) *SN ) *EXP < ALF*X ) + ! ( X ( 3 ) *ALF**2« '2.6 
2-2*0*K!4)*ALF#BET-K{3)#BET#*2<>0) #CS + ! K. ( 4 ) *ALF**2 .0+2 * 0*K. < 3 > *ALF#.F*B 
3ET-K(4)*BET**2.0)*SN)*EXP !-ALF*X) B 1 

UX=F66»<N0+D11*WXX/R-A12*W/R-F16#T/<F66*2*0*PI*R)-2*0«F16*D16*WXX)XX/ 
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1(R*F66))/D 

VX= < T*A1 l/( 2 #0*Pi *R )-F16*N0+WXX# £ 2 •0*D16*D-F16#D11#F66 + 2#0*D16*F1 , 6 
1**2#0 ) / (R*F66 ) +W*A12#F16/R ) /D 1 


Ed )=ux-z*wxx 

E £ 2 )= < 1 .0-Z/R+ £Z/R>##2«0 ) *W/R 
E ! 3)= ( 1.0+Z/R>#VX 
DO 2 1=1 »3 

STR £ I ) =0* 0 

DO 2 J= 1 9 3 

SIR £ I >=STR< I )+Q( I *J)#E<j ) 

2 CONTINUE 

WR ITE (3 *300 ) X*Z »NO*P»T »STR ( 1 ) 9STR< 2 ) *STR ( 3 ) »W»WXfWXX 
300 FORMAT (lH ♦F5»2»lX9F6«3»lX»F7 0 l*2X9F7*l»2X>F7.1t2X»Ella4«2X»Ell,4'» 
52X*£11*4»2X»E11«4*2X>E11«492X*E11#41 
X=FLOAT(N)»2*5/50»0 
1 CONTINUE 
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tie 


CALL FPLOT(-290«0»0«0) 

X = 0*0 

DO 31 N=1 ♦ 201 
CS-C05 ( BET#X I 
SN=SIN£BET*Xl 

W- < K ( 1 ) *CS+K < 2 ) *SN ) *EXP( ALF*X > + t K I 3 ) *CS+K < 4 > #SN ) #EXP ( -ALF*X ) +W0 
WX = U K( 1 >*ALF+K< 2 ) *BET )*CS+ ( K£ 2 ) *ALF~K.l 1 ) *BET > #SN > *EXP £ WLF*X ) + t(K{ 
84)*BET-K(3>#ALF>*CS-<K<3 ) #BET+K ( 4 ) *ALF ) *SN ) *EXP £-ALF*X ) ; 

WXX= ( < K ( 1 ) * ALF**2 •0+2a0*& ( 2 ) *ALF*SET-K< 1 ) *BET*#2 • 0 ) #CS + £ K ( 2 ) *ALF** 
12 #0-2 » 0#K ( 1 )*ALF*BET-K(2 ) *BET**2 # 0 ) *SN ) *EXP £ ALF*X ) + < { K £ 3 )*ALF##2*0 
2-2#0*K<4)*ALF*BET-KO>*BET**2.0) #CS + ( K < 4 J *ALF**2 « 0+2 »0#K { 3 J^ALF^B 
3ET-K( A)#BET**2.0)*SN)*EXP£-ALF#X) B :j 

UX=F6 6# £ N0+D1 1 #WXX/R-*A1 2 #W/R-F 16# T/ ( F66#2 «0#PI#R J -2#0#F16#D16*WX>i/ 
1£R#F66M/D .j 

Pi-CfTT^tS, VX=(T#A11/£2#0#PI *R )-Fl6*N0+WXX* £ 2 # 0#D16#D-F 16*D1 l#F66+2 *0#D 1 6*F16 

SCCTiod 1**2*0 ! / (R*F66)+W*A12#F16/R) /D 
E( 1>=UX-Z*WXX 

E(2) = ( 1#0“Z/R+(Z/R) ##2*0 ) #W/R 
E ( 3 >= ( 1 .O+Z/R )*VX 
DO 32 1=1*3 
5TR < I )=0,0 
DO 32 J = 1 » 3 

STR( I )=5TR( I >+Q( I » J)#£( J ) 

32 CONTINUE 

DO 73 1 = 193 

STR < I ) =STR £ I ) /1000*0 
73 CONTINUE 

CALL FPLOT<09X»STR£1) > 

X=FLOAT(N) #2.5/200*0 
31 CONTINUE 

jt CALL FPLOT£ l90*090»0> 


Z*Z- 0.02 
21 CONTINUE 
CALL EXIT 
END 
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The tube considered in this example was a circumferentially 
wrapped tube loaded to failure in simple axial compression. 

The expression used to calculate the ultimate failure stress 
was simply 

<Tult. = Pult. 

A 

and this gave a value f or <T X of -13*153 psi. The computer 
program was then used to calculate the stress distribution 
in the tube using the ultimate load as input and it gave the 
same results in the test section of the tube. It also showed 
that the tube had been reinforced sufficiently at the ends to 
overcome the large stresses that were induced there. As can 
be seen from the printout, the stress distribution across the 
thickness was also fairly uniform. 
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In this example, the tube was also wrapped circumferentially 
but it was under a more complex loading state. The tube was 
loaded in axial compression and internal pressure and the 
values of the load read into the computer program occurred at 
failure. The induced bending stresses were quite significant 
near the end of the tube that these rapidly converged to a 
uniform stress state. 
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APPENDIX 0 


THE EFFECTS OF TEMPERATURE ON THE MATERIAL PROPERTIES 
OF FIBREGLASS/EPOXY COMPOSITES 

B.l INTRODUCTION 

With the increased use in recent years of advanced com- 
posites by the aerospace industry for the design of primary 
load carrying structures, ever increasing consideration has 
been given to the effects that varying environmental con- 
ditions may have upon the integrity of these materials. Some 
of the most common of the advanced composites are those made 
using glass fibres embedded in a polymer matrix; an example 
of which is Scotchply 1002. Although some references exist 
on the effects of temperature on these materials it is de- 
sirable to have specific data on the specific material in 
question. This study has been performed in order to evaluate 
the effect of various temperatures on the elastic constants 
and burst strength of circular cylindrical tubes made using 
Scotchply 1002. The results of this investigation are re- 
ported in the following sections. 

B.2 Test Procedure 

All sixteen of the tubes used in this investigation were 
tested under internal pressure and subsequently loaded to 
failure. The various phases of the testing procedure are 
shown in Figures B-j {©'. 8J. Each of the eleven strain gauged 
tubes was subjected to two separate loading constraints; the 



first being free expansion under internal. pressure (ie*, both 
£ x and €y being non-zero) and the second being clamped ex- 
pansion under internal pressure such that £ x = 0 with, how- 
ever, €y remaining non-zero. By using the data obtained 
from the axial and circumferential strain gauges when the 
tubes were loaded as stated above, it was possible to solve 
the constitutive equations for the elastic constants of the 
material (ie., %22j ^ 12 and V 2 i). Section 3 describes 

in detail the method of solution used. Figures 8.4 through 
B.|0 show the pressure vs, strain curves from which the data 
was obtained to perform the calculations. 

In-situ test temperatures were varied over the range of 
-80*F to 250* F so that the effect of test temperature on the 
elastic constants of the material could be studied. Once 
this data was obtained the specimens were then loaded to 
failure under internal pressure in order to evaluate the 
effect of temperature on the transverse burst strength of 
the material. 

For these tests, temperatures in excess of ambient room 
temperature were obtained using a Blue "M" Oven, whereas 
temperatures below ambient were achieved using a Tenny Junior 
environmental chamber. In order to ensure that consistent 
in-situ test temperatures were obtained, each of the specimens 
was equ%)ed with thermocouples which were monitored at all 
times during each test. The temperature of each tube was 
allowed to stabilize at the desired value before testing be- 
gan. This nominally required 0,75 hours for each test. By 



using the Blue "M" Oven and the Tenny Junior Chamber it was 
possible to achieve a temperature stability of + 1*F during 
each test run. 

B.3 Basic Equations 

For any anisotropic material it is possible to write the 
constitutive or Hooke's Law relations in the following form: 


r* *1 
Ox 

k 

Q11 

$12 

Ql6 

<r y 

= 

Q12 

$22 

$26 

^xy 


Q16 

$26 

$66 


'-x 

£y 

€ xy 


( 01 ) 


where the superscript, k, denotes the k 1 -’^ lamina of the ma-r 
terial. Since all of the tubes used in this investigation 
were 3 ply, 90 ° laminates the superscript may, in this case, 
be omitted. 

The components of the stiffness matrix, [Q] are given 
by the relations: 

Qll = Qncos + 2(Q]_2 + 2 Q 55 ) sin 2 -«>cos 2 -@- + Q 22 sin 

Q 22 = Qnsin + 2(Qi 2 + 2Q66) sin 2 -©- cos 2 -e- + Q 22 cos 

Ql 2 = (Q 11 + Q 22 - 4Q66) sin 2 -©* cos 2 -©-+ Q 55 (sin^-o-+ cos4^) 

Q 56 = (Q 11 + Q 22 - 2Qi 2 - 2Q££) sin 2 -©-cos 2 -©- + Q 66 (sin^ + cos^-©) 

„ 1 

Ql6 = (Q11 - Ql2 - 2Q55) sin-©- cos 3 -©- + (Qi 2 - Q22 + SQ^sin^-©- cos •©■ 

Q 2 6 - (Qll “ Ql2 ~ 2Q 66) sin^©-cos 4 x+ (Q12 - Q22 + 2Q66>sin^- cos 3 

where the Qpj (i,j = 1 , 2 , 6 ) are the components of the stiff- 
ness matrix for the lamina coordinate system, and may be ex- 


( 02 ) 



pressed as: 


Qll = En/( 1 - V12V21) 

Q22 = E 22 /(l- Vi 2 V2l) 

Ql2 = v 21 e u/( 1 - V 12 v 2l) = V 12 E 22/( 1 “^12 V2l) (B 3 ) 
Q66 = G 12 


Ql6 = Q26 “ 0 


Thus, with a wrap angle of 90° as for the tubes being 
considered here, equations (02) reduce to: 

Qll = Q22 = 322/(1 “ ^12^21) 

§22 = Qll = E ii/(1 - V 12 V 2 i ) 

§12 = Ql 2 = V 21 En/(1 - VL2^2l) (B4) 

§66 = Q66 = G 12 
Ql6 = Q26 = G 

The constitutive equation (01) can therefore be expressed 
in matrix form as: 
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§22 Ql2 0 
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Q12 Q11 0 
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0 0 Q66 
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Equation (85) can now be expanded taking into consider- 
ation the two separate loading cases: 


<T X 1 * E 22 6xl +V12 E 22 fcyl (B6.a) 

3 3 

°yl = Vi 2 E 22 € X 1 + Ell € y l 

V~ T 


(86. b) 



and, 


^xyl - g 12 ^xyl 


( B6 . c ) 
( B6 . d ) 


(B6.e) 
CB6.f ) 


where the subscripts 1 and 2 denote the free expansion and 
clamped loading conditions ^respectively , and P = 1 “ V 12 ^21* 


By setting T y -j_ = 0“y2» equation (B6) 

can be reduced to 

the set: 

O’ xl = e 22 € xl + %2 E 22 Gyl 

(07. a) 

3 .0 


O’yl = ^12 e 22 G xl + E ll^yl 

(B7.b) 

3 3 

<Ty2 -Cyl = E 11 €y2 

(B7.c) 


3 

which may be readily solved for E22» ^12 and >>21 as 

follows, 

E H = <*>2 (B$) 

^ e y2 

E 22 _ jfel ' q'yl€yl (€ y 2 ~ £yl) 

0 ^xl \ ^ xl £y2 



°x2 = ^12 E22 £ y2 

r~ 

<r y2 = e x1 e y2 

0 

°*xy2 “ g 12 ^ xy2 


V 12 = < Tyl (€ y 2 ~ €yl) €xl 

^xl £ xl &y2 -^“yl€ y i(£y2 - 6yl) 


( 010 ) 



V 21 = (^y2 - €yj) - - . .. (Bll) 

£ xl 

It should be noted that the condition thatO*yi = 0y2 is 
easily accomplished experimentally by ensuring that under 
both of the loading conditions the tubes are loaded to the 
same internal pressure. 

The above calculations were performed for each tube and 

the results are listed in Table IV. The geometric properties of each of the 
tubes tested can be found in Table B.l. 

B«A Discussion of Experimental Results 

As can be seen from the graphs in Figures B*ll to S-^the 
testing temperature has little effect on the value of E^]_, 

The calculated values are well within + 10$ of the average 
value of about 5*9 x lO^psi. It was expected that the range 
of temperatures investigated would not affect the properties 
of the brittle fibres as much as it would affect those of the 
ductile epoxy matrix. Therfore, since is primarily a 
function of the fibre modulus it should not change very much. 

The test temperature had a greater effect on the value 
of E 22 since it is primarily a function of the modulus of the 
matrix. At high temperatures the matrix is very plastic and 
strains easily so that E 22 takes on a very low value. At 
lower temperatures the matrix becomes brittle and strains 
less than it would at higher temperatures. This increased 
slope of the stress-strain curve results in a higher value 
for E 22 . 



The same reasoning applies to the effect of temperature 
on the transverse burst strength of the material. Since it 
strains easier at higher temperatures, fracture will occur 
at a lower value of internal pressure and the opposite is 
true for the lower temperatures. 

B.5 Conclusions 

The effects of testing temperature on the transverse 
burst strength and the material properties of a fibreglass/ 
epoxy composite have been reported in this /ippendi'x. It has 
been shown that variations in the temperature over the range 
being investigated have pronounced effects on the properties 
of the matrix material but do not affect the fibres signifi- 
cantly. Since it is the matrix material which is most af- 
fected by the temperature changes, the effects on other mech- 
anical properties such as shear strength and modulus which 
are directly related to the matrix should be investigated 
fully in order to obtain a more complete understanding of 
the changes taking place in the composite at different 
temperatures. 
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FIGURE BJ 

INSTRUMENTATION AND A TUBE READY FOR TESTING 
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FIGURE B.2 


CLAMPING APPARATUS AND A FRACTURED TUBE 




FIG. B3 CHAMBERS FOR THE COLD AND HOT 
TEMPERATURE TESTS 


STRESS (RS.L) 
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APPENDIX C 


THE EFFECTS OF LENGTH OF POST CURE ON THE MATERIAL 
PROPERTIES OF FIBREGLASS/EPOXT COMPOSITES. 

C.l INTRODUCTION 

Many fibreglass/epoxy components can be manufactured 
from preimpregnated (prepreg) sheets of oriented fibre im- 
pregnated with un cured epoxy resin. Data sheets provided 
by the manufacturer of the prepreg tapes usually give a 
specified length of time as well as the temperature and 
pressure at which the component must be cured to achieve 
optimum properties. However, these specifications rarely 
give the effects of a deviation from the normal cure time 
on the properties of the composite. 

Here, the results of a preliminary investigation on the 
variation in certain mechanical properties which result from 
variations in the length of post cure time will be reported. 

C . 2 Test Procedure 

All of the tubes were tested under internal pressure. 

Two tubes from each batch of three were equipped with an axial 
and a circumferential strain gauge. The elastic constants 
(E^q, E22»Vi2> andV2l) of the material can be solved for 
when the responses from these strain gauges to two different 
loadings are substituted into the constitutive equations (see 
Appendix 13 for details). The two types of loadings 

are free expansion under -internal pressure (6 X ^0 and€yj*0) 



and clamped expansion under internal pressure where the tube 

V 

is clamped in such a way that € x = 0 but fcy $ 0. 

The required pressure vs. strain curves were obtained 
for each strain gauged specimen (Fig.C.I to C.4 ) and then the 
tubes were burst in order to find the effect of length of post 
cure on the ultimate transverse burst strength of the material 
(Fig.C.5). 

A general view of the testing apparatus is given in 
Figures B. I andBc2. 

All of the tubes used in these tests were wrapped with a 
fibre angle of 90°, ie. , the fibres were aligned in the cir- 
cumferential direction. 

The equations used in the calculations are the same as 
those derived in the previous Appendix. 

C.3 Conclusions 

The various elastic constants and ultimate transverse 
strength were calculated from the strain gauge responses for 
the two different loadings. These results are shown in 
Table V and are plotted in Figures C* 6 +0 Ch9. 

As can be seen from the graphs, the length of post curing 
time has little effect on the values of E]_p The calculated 
average values vary little from about 5.6 x 10^ psi. This 
was expected because is essentially the value obtained by 

multiplying the fibre modulus by the fibre volume fraction. 



The effect of the matrix modulus on Ell is negligible. Also, 
the lengths of post curing times considered in this investi- 
gation shouldn't have had any degrading or harmful effects on 
the fibres. For the 2 4 hour post cure case, there was quite 
a variation in the values calculated for En. This was de- 
termined to be a Poisson's ratio effect. The Poisson’s ratio, 
Vi2> is directly proportional to the difference in the strain 
Gy2 andGyX (a very small number) and is therefore very diffi- 
cult to determine accurately. The difference in the two 
values is approximately 0.2 which gives a + 33$ inaccuracy in 
the value of Vi2» 'Ell was directly proportional to Vi2 
(Ey]_ ='^12 E 22 /V 21 ) and this is the reason for the great vari- 
ation in 24 hours postcure. The quantity 3 - (1 - V 12 V 21 ) 

was a much more accurate number to determine and so the quanti- 
ties Eii/jj and E 22/3 were calculated so that the variation 
in the moduli would be reduced. 

The difference in post curing time has a more pronounced 
effect on the value of E22« As the length of post curing 
time increases, E 22 steadily decreases. This is also to be 
expected because the length of post cure affects the matrix 
more than the fibres and E 22 is more dependent on the modulus 
of the matrix. 

Using the same reasoning as above, the ultimate trans- 
verse strength, 5*2 u> increases a good deal as the length of 
post cure increases. 

From these two facts it is deduced that an increase in 



the length of post cure causes the matrix to become more 

ductile. 

The amount of post cure will probably affect other me- 
chanical properties of the composite more than was noted in 
this investigation. Some of these other properties would be 
the interlaminar shear strength and shear modulus, the trans- 
verse compressive strength and the fatigue strength. The 
effects o.n these properties will probably be quite pronounced 
since they are quite dependent on the state of the matrix. 
Other investigations should provide the effects of differences 
in post cure temperature as well as time in order to determine 
the optimum curing time for any given curing temperature to 
obtain the optimum value for any particular mechanical pro- 
perty of the composite. Curing conditions will probably 
vary with the property that is required and so, different 
components can be fabricated to optimize the specific me- 
chanical property that is required. 



TABLE C.l 



GEOMETRY 

OF TUBES USED IN POST CURE 

STUDY 

TUBE 

DESIGNATION 

R 

(IN) 

T 

(IN) 

11c 

90° 
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1.016 

.0295 

12c 

90° 

1.016 

.0295 

13c 

90° 

1.016 

.0295 

18 a 

90° 

1.016 

.0292 

18b 

90° 

1.016 

.0297 

18c 

90° 

1.016 

.0300 

20a 

90° 

1.016 

.0299 

20b 

VD 

O 

0 

1.016 

.0298 

20c 

0 

o 

i 

1.016 

.0292 

21a 

90° 

1.016 

.0294 

21b 

90° 

.1.016 

.0303 

21c 

90° 

1.016 

.0303 
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APPENDIX D 


COMPUTER PROGRAMMES FOR EVALUATING CUBIC TERMS M) FAILURE 
; PRESSURES FOR LAMINATED TUBES UNDER INTERNAL PRESSURE- 


WITH AND 'WITHOUT TORQUE 





EVALUATION OF CUBIC COEFFICIENTS 


IMPLICIT REAL*8 < A-H>0~Z ) ' ~ ' 

COMMON /AAA/Ell »E22,V12,V21 »G12 

DIMENSION PRRC4) ,ENXY<4) » DT I ( 4 I ♦ TR ( 3 1 3 ) S I G < A , 3 ) ,A < 4 » A ) ,B ( 4 ) * 

1 I R C 4 ) * I C C 4 ) 

* READ ( 5 ♦ 2 ) Ell »E2 2 *V12*V2 1 *612 

2 FORMAT (5D16. 4) 

READ (5,5) F1»F2»F11»F22*F12>F66 
5 FORMAT (6D13.7) 

. READ (5 > 10 ) . ( PRR.C I L»I=1.*.4J 

10 FORMAT (4D20 .5 ) ’ “ " 

READ (5*10) (ENXY ( I ) *1 = 1 * 4) 

READ (5,10) (DTI < I ) , I = 1 > 4 ) 

READ C 5 1 1 > LR 
1 FORMAT (110) 

TH= • 04000 ; 

RAD* 1*0200 
NL=4 

DO 40 1*1*4 - ' 

DT = DT 1(1) 

CALL CMAT(TR,DT»TH.NL»LR) 

SIGU » 1 ) =TR (1 ,1 ) *PRR( I ) / 2.D0 + TR ( 1 »2.) *PRR< 1 )+JR< 1 ,3 )*ENXYt I 

SIG( I *2 ) *TR ( 2 » 1 ) *PRR( I ! / 2 « DO+T R ( 2 » 2 ) *PRR < 1 ) +TR l 2 »3 ) * ENXY l 1 ) 

SIG( I t 3 ) =TR ( 3 » 1 J *PRR ( I } /2.D0+TR ( 3,2 ) *PRR< I ) + TR( 3,3) *ENXY( I ) 

.40 CONTINUE 

DO 50 1-1,4 

A ( I * 1 ) =3 • D0*S I G ( I » 1 ) *S I G ( I * 1 ) * S I G ( I »2) 

A C I . 2 >=3.D0*SIG( 1,1) *S I G C I , 2 ) *S I G l I * 2 ) 

AU*3)=3.D0*SIG( I»2)#SIG( I ,3)*SIG( I »3> “ ' 

AC I .4 >=3.D0*SIG ( I ,1 >*SI G ( I ,3 >*SIG< I ,3 ) 

B LI )=1.D0-CF1»SIG(I,1)+F2*SIG(I,2)+2.D0»F12*5IGCI,1)*SIG(1.2} 

1 +F1 1 *S I G ( I »1)*SIG(1»1)+F22*SIGC I » 2 ) *S I G ( I *2 )+F66*SIG( I , 3 ) *S 1 G ( I , 3 

1 ) ) 

.50 CONT I NUE • _ 

CALL LNEQND C A ♦ 4 ♦ 4 ♦ 4 V B ♦ I R » I C » f E R T 

WRITE (6,60) 

60 FORMAT C 1 1 ' ,T10, 'F112 1 »T30, 'F221 1 ,T50, «F 266 * » T 70 , » F 1 66 ' , // ) 

WR I T E ( 6 * 70 ) <B(I)*I»1»4> 

70 FORMAT ( 4D20 • 7 J 

« STOP 

END ~~ ' " " ” 

•4765D+07 • 12 1 1D+ 0 7 .3363D+00 .0855D+00 .4230D+06 

~ * 307558 90-0 2 + , 23 440 1 00+00+ , 93 9 771 1D-U4+, 2269570D-01-, 638730 C)p-0 3+ « ? Ifr 74 1 flp-p y 

•67422D+00 ,674220+00 *601800+00 .6Q180D+0Q 

•OOOOOD+OO ,076500+00 ,076500+00 .OOOOOD+UQ 

.600000+02 , 60000D+.02 ;.,45000D + 02 .450000+02 



CALCULATION OF FAILURE PRESSURES 


SPRINTOFF 

IMPLICIT REAL *8 ( A-H >0-Z ) 

COMMON /AAA 1 FI 1 . F 2? .VI 2. V 2 1 .fil 2 . ...... 

? 

DIMENSION TR ( 3 »3 ) 

READ ( 5 » 2 > Ell »E22*V12»V21»G12 
FORMAT ( 5016.4 > 


5 

READ ( 5 » 5 ) Fl«F2»FlltF22*F12»F66 
FORMAT ( 6D 1 3 • 7 ) 

I R = 2 _ .. . 


ENXY = • 0765D0 
DT = 8 0 • DO 

Fll 2 = . 1 6670670-03 

F221=-.l057499D-02 
F266=-. 1372627D-02 

F] 66 = * 399 8 A 9 AD-0 3 , ....... • ■ - 

TH= • 040DO 
DP= • 005DO 
RAO-1. 0200 


AL I M= • 8D0 
NL = 4 

WP I TF t A , 1 n > DT .FNVY • 


10 

FORMAT { ' 1* » T 1 0 ♦ ' THETA= ' » F10.2»T30« 'NXY=' tF12*4*//) 

CALL CMAT ( TR»DT »TH»NL.LR ) 

P R R = 0 • DJ2 


40 

CONTINUE 

SIG1=TR( 1«1 )*PRR/2.D0+TR< 1 ,2 )*PRR+TR(1 ,3)*ENXY 

SIG2=TR( 2 .1 l*PRR/2« 00+T R (2 *2 )#PRR + TR(2»3)*ENXY „ ... 



SI 66 = TR ( 3 1 1 ) *PRR /2 .DO+T R ( 3 * 2 ) *PRR+TR ( 3 ♦ 3 ) *ENXY 

VALF=F1*SIG1+F2*SIG2+F1 1 *S I G 1*S I G 1 +F 2 2*S I G2*S I G2 +F66 *5 I 06*S l G6+ 

1 2.D0*F12*SIG1*SIG2+3.D0*F112*S I61*SIG1*SIG2+3.D0*F221*SI62*SI62* ... 



1 SIG1+3.D0*F266*SIG6*SIG6*SIG2-1 .D0+3.D0*F166*S I Gl*S I G6*S I G6 
PR=PRR 

WR I T E ( 6 * 850 ) PR * VALF »S I G 1 ♦ S IG2 ♦ S I G6 


8 50 

FORMAT ( • • , *PR=' ♦ F 10 • 5 t 2 X ♦ * VALF = 1 » D 1 5 . 5 » 2X ♦ • S I G1 = ' » D1 5 . 5 ♦ 2 X ♦ * S 1 62 

1= » *D15.5 »2X» 'SI 66= ' *D15. 5 ) 

PRR= PR R + DP ... -- 


V 

100 

iff- 

c 

I F ( PR .GE. ALIM) GO TO 100 
GO TO 40 

CONTINUE * 


STOP 

END 



SUBROUTINE CM AT (TR»DT ♦ T H * NL » LR ) 

IMPLICIT REAL* 8 { A-H ♦0-Z ) 

COMMON /AAA/E11 *E22*V12»V21*G12 

DIMENSION T (4 ) «DTA<4 ) »T A ( 4 ) ,TR ( 3 *3 > *TRI < 3 * 3 ) *Q( 3 »3 ) »H ( 5 ) ♦ 

1 OQ(3»3»4) *A(3»3) »AA(3»3) ♦ LL < 3 ) » MM ( 3 ) * QR < 3 ♦ 3 ) . W < 3 ) 

DO 88 K = 1 *NL ’ _ - 

- 



n n 


00 ( 1 ♦ 2 ♦ K ) = 014*S2*C2+Q( 1 *2 MMS4+C4) 

00 { 2 ♦ 1 ♦ K } = 0Q ( 1 * 2 ♦ K > 

— QO(3»3»K) = Q22*S2*C2+Q( 3 » 3 > * ( S4+C4 ). ... 

QO ( 1 ♦ 3 ♦ K ) = (012 +C2+Q23* S2 ) *S 1+C 1 
Q0(3*1*K) = 00(1 »3*K) 

QQ t 2 ♦ 3 » K ) = tQ1?»S2+Q23*C2)«Sl»C1 

9 QO ( 3 * 2 ♦ K ) = QQ ( 2 * 3 *K ) 
v DO 10 J = 1 » 3 

, DO. 10 1=1*3 ... „ 

A ( I » J ) =0 • DO 
DO 10 K=1 ♦ NL 

1 A_LL t J ) =QQ ( I »J*K)*(H(K+1 ) — H ( K ) > ± A(I*J) 

10 AA ( I* J) = A( I * J) 

SUBROUTINE MINV ENTER AA AND REPLACE IT WITH THE INVERSE MATRIX 

MM,LL ARE WORK VECTORS FOR SUBROUTINE M I NV„ J MATH I X INVERSION). _ 

CALL MINVRD(AA#3*3,DET*IER*LL»MM) 

DO 12 1=1*3 

DO 12 J = 1 *3 __ __ 

12 QR ( I * J) =QQ( I » J.LR ) 

CALL MAMP 1 D ( T R » 3 » 3 »QR * 3 » 3 » 3 * 3 » 3 * W * 3 ) 

1 CALL MAMP1D(T.R»3 .3_t.AA»3 t 3»3»3»,3 »W*3 ) .„ _ 

RETURN J ----- 

END 

SDATA 

SDATA _ ' 

•4765D+07 *12110+07 *33630+00 .08550+00 *42300+06 

- . 307 65 8 ? D - 02 + . 2 34401 00+00+ • 9 3977 1 .1 D - 0 4+ . 2 2 6 9 5 7 0 D- 0 1 -*6387 3 0.0 Q - 0 3 + . 2 142418 D-02 


STATISTICS ~ 136 CARDS READ*- 150 LINES PRINTED 


